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INTRODUCTION 


Scope 


This  thesis  discusses  a  special  purpose  digital  computer 
for  the  arctangent  function.   Reasons  for  this  choice  will  be 
explained,  and  the  computer  design  "will  be  based  on  an  upper 
polygonal  approximation  with  uniform  Chebyshev  norm  for  the 
error. 

A    special  purpose  digital  computer  (SPC)  executes  a  small, 
special  set  of  instructions  and  solves  only  a  restricted  number 
of  problems,  perhaps  one  problem.   On  the  other  hand,  a  general 
purpose  digital  computer  (GPC)  solves  a  large  set  of  problems. 
In  past  years,  both  types  of  digital  computers  have  been  con- 
structed.  The  GPC ' s  are  almost  exclusively  employed  in  computing 
centers  of  industry,  universities,  and  business  corporations. 
However,  SPC ' s  came  into  existence  long  before  GPC's.   Familiar 
examples  of  SPC ' s  are: 

1.  The  switching  system  in  an  automatic  telephone  exchange. 

2.  Jacquard's  punched  card  control  for  a  loom. 

3.  Numerical  control  of  drilling  and  milling  machines. 
SPC ' s  are  especially  important  in  guidance  systems  of  missiles 
and  satellites.   Economical  and  technical  considerations  demand 
SPC • s  for  many  applications. 


The  Computer  as  an  Abstract  Notion 

The  term,  "computer",  does  not  connote  the  metal  parts,  com- 
ponents, circuits,  functional  units,  and  so  on.   The  underlying 
concept  of  a  digital  computer  is  the  proven  algorithm.   The  dis- 
tinction between  a  GPC  and  an  SPC  may  be  narrowed  by  the  follow- 
ing observation:   A  general  purpose  computer  is  a  collection  of 
special  purpose  computers  (Pigs.  1  and  2). 

Nonerasable  and  Erasable  Memories 

There  are  two  types  of  memories: 

1.  The  nonerasable  memory  has  memory  cells  whose  contents 
remain  unchanged  during  the  computation  process.   There 
is  no  read-in  of  new  information,  and  the  read-out  is 
nondestructive.   Everyday  examples  are:  dictionary, 
photography  album,  phonograph  record,  papercard  with 
punched  holes.   Nonerasable  memories  are  likewise  de- 
noted as  permanent  memories  or  as  read-only  memories. 

About  ten  years  ago,  the  capacitive  and  the  magnetic 
memory  types  aroused  interest.   Their  features  and 
operation  are  described  in  (2,  3,  !).).   A  general  survey 
of  nonerasable  memories  is  given  in  (1).   They  are  used 
to  store  subroutines  and  other  kinds  of  information. 

2.  The  erasable  memory  has  memory  cells  with  the  property 
that  their  contents  can  be  replaced  any  time  by  new 
entrie  s. 
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Fig. I..  Symbol   for  an    SPC. 
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Fig.  2.  GPC   as    a  collection    of   interlinked    SPC'S. 


Figures  3  and  if  show  schematically  these  two  types  of 
memory  cells. 

AN  ARCTANGENT  SPECIAL  PURPOSE  COMPUTER 

Problem 

Our  objective  is  to  design  an  SPC  that  calculates  bhe  arc- 
tangent function  with  a  preassigned  accuracy.   After  a  brief 
discussion  of  the  power  series  method,  some  properties  of  the 
arctangent  function  will  be  indicated.   As  it  is  a  transcendental 
function,  the  value  for  a  certain  argument  is  obtained  in  an 
infinite  number  of  steps  by  means  of  power  series  such  as: 

1,1,-  ? 

arctan(x)  =  x  -  —  xJ  +  —  yp    -  .  .  .  ,  x  <  1       (1) 

3     5 

11  i   x  -i   1  z 

arctan(x)  =--xx+-xs>--   x~p  +  .  .  . ,  x  >  1   (2) 

2         3       5 
Equation  (1)  ie  derived  in  (6)  by  means  of  Taylor  series.   Equa- 
tion (2)  is  obtained  from  equation  (l).   This  is  done  in  the  in- 
troduction of  reference  (7).   Reference  (2)  contains  a  12-digit 
tabulation  of  the  arctangent  function.   The  limitation  to  exe- 
cute only  a  finite  number  of  steps  gives  rise  to  questions  of 
approximation,  speed  of  convergence,  and  error  measure. 
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Fig. 3.  Non-erasable  memory  cell. 
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Fig.  4.  Erasable  memory  cell. 


Convexity,  Symmetry,  and  Monotonicity 

Let  f(x)  be  a  function  defined  in  the  open  interval 
I  =  (a,  b)  ,  and  let  I-,  =  (x-,,  Xp)  be  an  open  interval  such  that 
I]_  cl   I.   g(x)  is  a  line  through  the  points  jx-j,  f(x-,)j  and 
[x-2*    ^(x2^1  *   f^x)  ^-s  convex  in  I,  if  for  all  x  €  1-^  holds 
fix)    -  g(x)  <  0,  whichever  1^  is  chosen  (Fig.  5) •   Dually,  f(x) 
is  concave  in  I,  if  for  all  x  €  1-^  holds  f(x)  -  g(x)  ">   0, 
whichever  I-]_  is  chosen  (Pig.  6)  . 

Convex  functions  are  also  called  sublinear  functions,  and 
concave  functions  are  called  superlinear  functions .   If  a  func- 
tion is  twice  differentiable,  convexity  corresponds  to  positive 
curvature  and  concavity  to  negative  curvature.   A  short  and  ele- 
gant theory  on  convex  functions  is  given  in  Emil  Artin's  mono- 
graph (5)  . 

Analyticity  and  concavity  of  the  arctangent  function  in  the 
open  interval  (0,  oO)    imply  that  it  has  negative  curvature 
throughout  (0,  oO) ,   By  Rolle ' s  Theorem,  the  minimum  of  curvature 
lies  in  (0,  oO ) .   The  calculation  of  Appendix  A  shows  that  the 
minimum  occurs  at  x  =  0.831. 

The  arctangent  function  is  skewsymmetric  with  respect  to 
the  origin;  it  is  an  odd  function.   Therefore  only  the  open  in- 
terval (0,  cxo)    is  relevant  to  us. 

With  the  argument  x  increasing  throughout  the  open  interval 
(-00,  00) ,  the  arctangent  function  is  increasing  too,  and  there- 
fore monotonic  in  (-00,  00  ) .      The  arctangent  function,  of  course, 
is  considered  as  single  valued,  defined  in  the  value  range 


Fig.  5.  A  convex  function    f(x) 


Fig.  6.  A   concave  function    f(x). 
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Error  Measures 

Different  ways  to  measure  the  quality  of  approximation  are: 

(i)  error  area   /  |e(x)|  dx 
X 

(ii)  least  square  area  J     e(x)    dx 
X 

(iii)  maximum  of  absolute  error  value,  max  |e(x) 

X   '     ' 

where  e(x)  is  the  error  function.   Error  criteria  (i)  and  (ii) 

have  a  great  disadvantage.   They  are  integral  error  measures  and 

do  not  take  account  of  local  errors.   Gibbs  phenomenon  in  Fourier 

series  is  a  case  in  point  of  this  disadvantage.   Error  criterion 

(iii)  notes  local  errors.   They  are  bounded  by  the  preassigned 

maximal  error  ©max«   This  error  measure,  the  Chebyshev-norm, 

will  be  employed. 


POLYGONAL  APPROXIMATION 


Definition 


A  way  to  approximate  a  continuous  function  f(x)  (Fig.  7) 
is  to  use  an  approximating  polygon  (Pq  P-^  P2  .  .  •  %  pn+1^  * 
The  vertices  of  this  polygon  coincide  with  points  of  the  func- 
tion such  that 

P0  =  {xQ,  f(xQ)}  ,  ?1   =    {xlt    f(Xl)},  .  .  ., 

PN  =    1XN'  f^xNM  '       PN+1  =  ixN+l'  f^xN+l^l  ' 

The  purpose  of  this  thesis  is  to  approximate  f(x)  = 
arctan(x)  in  the  interval  (0,«o)  by  a  polygon  (Pig.  8).   As  this 
function  is  concave,  the  approximating  polygon  is  also  concave. 

The  arctangent  function  has  two  inflection  points,  one  at 
zero,  coincident  with  Pq,  and  one  at  infinity,  coincident  with 
P^o  .   Hence  the  polygon  is  denoted  by  (Pq  P-,  P«  .  .  .  P-vr 
]?N+1  P^  )  .   In  functional  notation,  the  polygon  (Pq  P-,  P2  .  .  . 
P^u  pn+i  Poo  )  is  denoted  by  g(x).   g(x)  is  a  continuous,  monotonic, 
and  single-valued  function,  defined  in  the  domain  0  -£■   x  ^  ©o 
with  the  range  0  ^  y  <  +  ti/2  . 

The  orthogonal  projection  of  polygon  vertices  Pq,  P-j_,  .  .  ., 
PN>  PN+1>  poo  yields  the  end  points  AQ,  A-j_,  .  .  .,  AN,  AN+1,  kOQ 
of  the  intervals  IQ,  I,,  .  .  .  IN,  %+]_*  I©o  on  tiie  x-axis  such 
that  I0  =  (A0  Aj),    I1   =    (A1  A2),  .  .  .,  IN  =  (AN  AN+1)  and 
loo  =  (%+1  A^)  (Fig.  8). 

In  order  to  calculate  the  value  for  the  argument  x,  the 
computer  has  first  to  decide  in  which  interval  In  the  argument 
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Fig.  7.  Polygonal  approximation    g(x) 
of  a  continuous   function    f(x) . 
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Fig. 8.  A  polygonal  approximation  of  the 
arctangent   function. 
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x  lies.   If  x  €  I  ,  the  value  at  x  is 
x  •  mn  +  bn  =  y 
where  mn,  bn  are  stored  as  nonerasable  information  in  the  com- 
puter memory. 

The  following  subchapters  treat  in  detail  the  generation 
of  various  polygonal  approximations  for  the  arctangent  function. 
All  of  them  are  upper  polygonal  approximations,  i.e.,  polygonal 
approximations  that  satisfy  the  inequality  g(x)  -  f(x)  ^  0.  As 
the  arctangent  function  f(x)  =  arctan(x)  is  concave  in  interval 
(0,  ©o  ),  these  upper  polygonal  approximations  consist  exclusively 
of  tangent  lines. 

Tangent  Line  Intersection  Method 

This  method  begins  with  two  tangent  lines  of  arctan(x)  with 
points  of  tangency  at  x  =  x^q  =  0  and  x  =  x^00  =  oo .   Subscript 
"t"  denotes  the  tangent  line  of  arctan(x)  at  the  corresponding 
arguments  x.  These  two  lines  intersect  in  the  point  ]_PGo  (Fig* 
10)  with  the  argument  ]_xeo*   Subscript  "1"  denotes  the  first  ap- 
proximation stage,  and  "e"  tells  that  at  argument  x,  the  maximum 
error  of  this  polygon  approximation  occurs. 

This  first  approximation  divides  (0,  oo )  in  two  intervals, 
jIq    and  ]_Ioo  >    adjoining  at  the  dichotomy  point  ixe0'      ^e    ^an" 
gent  lines  of  arctan(x)  at  x  =  0  and  x  =  oo  are  called  half 
tangent  lines,  since  only  one  side  of  their  point  of  tangency 
is  considered;  if  both  sides  are  considered,  it  is  called  a  full 
tangent  line .   The  number  of  full  tangent  lines  employed  in  a 
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polygonal  arctangent  approximation  is  signified  by  N  and  the 
total  number  of  tangent  lines  by  M.   This  terminology  holds 
throughout  this  thesis. 

Two  consecutive  tangent  lines  of  arctan(x)  of  the  k 
approximation  stage  with  arguments  ^xt  •  and  ^xt  -+1  are  called 
adjacent.   Hence  tangent  lines  for  x  =  xtQ  and  x  =  xt  00   are 
adjacent  in  the  first  approximation  stage,  but  not  in  higher 
stages. 

It  was  discussed  how  to  arrive  at  the  first  approximation 
stage.   In  order  to  obtain  from  the  k^   stage  the  (k+1)    stage, 
the  points 

(kxeO>  kyeoj  >    ••*>    {kxen>  kJen]    >       •>    ^kxeN(k)>  kyeN(k)} 
are  constructed  by  intersecting  adjacent  tangent  lines  of  the 
k^*1  stage.   The  maximal  error  arguments 

kxeO'  kxel'  "•'  kxen'  •"'  kxeN(k) 


keO   k+l^tl'  k  el   k+l*t3'  ""'  ken    ( k+1)  e( 2n+l) ' 

(1) 
Now  the  tangent  lines  of  arctan(x)  for  the  arguments 

k+lXtl'  k+lXt3'  '"'    k+lXe(2n+l)' 
are  constructed. 

The  tangent  lines  of  the  k.^  stage  become  tangent  lines  of 

the  (k+l)th  stage  such  that 

kxtO  =  k+lxtO'  kxtl  =  k+lxt2'  •*'  kxtn  =  (k+l)xe(2n)' 

(2) 
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Equations  (l)  and  (2)  form  the  set  of  tangent  arguments  of  the 
points  of  tangency  for  the  (k+1)    stage.   The  relations  between 
tangent  arguments  and  maximal  error  arguments  are  represented 
in  Pig.  9. 

Since  each  maximal  error  argument  i,x   generates  two  new 
error  arguments  j£+ixe(2n)  anc*  k+lxe(2n+l)  ^-n  t-'ie  nex"t  stage, 
the  number  of  maximal  error  arguments  of  the  kth  stage  is 
NQ  =  2k~1    . 

(For  the  first  stage  holds  N   =  1.)   As  tangent  arguments 
and  error  arguments  follow  alternately,  there  are 

Nt  =  2k~1   +    1 
tangent  arguments  in  the  k    stage.   After  having  dropped  initial 
half  tangent  line  at  t^x.q  =  0  and  ended  half  tangent  line  at 
kxtoo  =  °'  "there  are 

N  =  N(k)  =  Nt  -  2  =  2k_1  -  1 
full  tangent  lines  for  the  k""1  stage.   This  shows  that  the  number 
of  tangents  employed  in  this  type  of  polygonal  approximations 
grows  exponentially  with  the  number  of  stages. 

An  important  problem  of  polygonal  approximation  is  to  re- 
duce the  number  of  tangent  lines  for  a  required  accuracy  in 
order  to  get  a  fast  computing  SPC  and  to  economize  on  memory. 
This  approximation  method,  however,  does  not  exhaust  suffi- 
ciently the  range  of  possible  polygonal  approximations.   Also, 
the  maximal  errors  are  not  uniform.   The  only  advantage  of  this 
approximation  is  that  it  is  easy  to  calculate. 


Ill 


E      Q. 
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Quasi-equal  Error  Method 

The  initial  stage  of  approximation  is  the  same  as  the  pre- 
ceding one  (Pig.  11) . 

The  second  stage  is  obtained  by  drawing  a  tangent  line  2^1 
that  contacts  the  arctangent  curve  at  the  point  {2xtl'  2^'tll'  * 
2x  ,  lies  in  the  interval  (0,  00   ).   This  tangent  intersects  the 
two  tangent  lines  of  the  initial  stage  in  the  points  {2xe0>  2^"e0 } 
and  {2xel>  2xel\  *   Therefore  the  approximating  part  of  the  tan- 
gent line  2^1  i-s  determined  by  the  interval  (2xe0'  2xel^  *  ^e 
argument  2xtl  ^s  c^-osen  such  that  the  maximal  errors  at  the  end- 
points  of  interval  (oxe0f    2xel^  ar>e  e1ua1'   Thus  the  second  stage 
yields  an  approximation  with  uniform  maximal  errors,  i.e.,  in 
the  uniform  Chebyshev  norm. 

This  procedure,  when  applied  to  all  pairs  of  adjacent  tan- 
gent lines  of  the  second  stage,  leads  to  the  third  stage.   The 
entire  procedure  of  the  second  stage  can  be  repeated  arbitrarily. 
It  generates  new  intervals  that  are  always  shorter  than  the 
preceding  ones.   As  numerical  computation  demonstrates,  the 
stages  of  higher  order  than  two  do  not  yield  uniform  maximum 
errors  at  the  endpoints  of  all  intervals. 

The  number  of  tangent  lines  for  the  k""1  stage  is 

k(k+l) 

Nt(k)  =  +  2 

t        2 

Nj-(k)  grows  exponentially  as  k  increases.   Also,  the  calculation 

to  these  approximations  is  not  direct;  it  involves  iteration. 
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2At0         2*eO        2*t| 


Fig.  10.  The  tangent  line  intersection  method. 


2AtO         2AQ0  2Atl 


Fig.ll.    The   quasi -equal    error  method 
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The  "Quasi-equal  Error  Method"  yields  by  transition  to  higher 
stages  smaller  and  smaller  errors.  In  the  following  chapter, 
the  "Polygonal  Approximation  with  Uniform  Chebyshev  Norm"  has 
not  only  this  quality;  it  has  uniform  maximum  error  at  the  end- 
points  of  the  intervals  IQ,  1-^,  .  .  .,  IN,  %+]_>  ^  oo  •  In 
addition,  polygonal  approximations  are  for  any  number  of  tan- 
gent lines  available. 

It  should  be  noted  that  the  polygonal  approximations  ob- 
tained by  using  either  the  tangent  line  intersection  method  or 
quasi-equal  error  method  are  all  upper  polygonal  approximations 


POLYGONAL  APPROXIMATION  WITH  UNIFORM 
CHEBYSHEV  NORM 


Definition 

If  an  approximation  is  subject  to  the  uniform  Chebyshev 
norm,  all  relative  maxima  of   e(x)   have  the  same  magnitude 
(Pigs.  12  and  13).   This  requirement  automatically  reduces  the 
number  of  tangents  used  for  polygonal  approximation.   The  third 
attempt  to  approximate  a  polygonal  approximation  can  be  de- 
scribed as  follows. 

Design  a  special  purpose  digital  computer  that  calculates 
the  arctangent  function  with  a  preassigned  accuracy,  using  an 
upper  polygonal  approximation  with  uniform  Chebyshev  norm. 
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Fig.  12.    Error    distribution   e(x)  of  a  polygonal   approximation 
with  uniform    Chebyshev  norm  . 


/My) 


1.57... 


Fig.  13.    Polygonal    approximation  with   uniform 
Chebyshev  norm. 
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The  Ch.ebysh.ev  Polygonal  Approximation 
with  Zero  Boundary  Error 


Zero  boundary  error  means  that  the  error  of  this  approxima- 
tion is  zero  at  the  inflection  points  of  the  arctangent  function. 
The  procedure  is: 

(i)  The  half  tangent  line  tg  of  the  arctangent  graph  at 
x  =  0  is  used  as  initial  approximating  line.   There- 
fore the  error  at  point  Tq  =  (x,yl  =  {°>°}  is  zero. 
From  here  the  point   x,y   travels  on  tangent  line  up- 
wards to  the  right.   The  error  e(x)  =  x  -  arctan(x) 
is  measured  simultaneously.   e(x)  increases  monotoni- 
cally  with  x.   Finally,  point  (x,y)  reaches  point 
P-]_  =  (x  -,  ,  Yei\   with  argument  x  =  x  ^.   At  this  posi- 
tion, e(x)  =  e(x  ■,)  is  equal  to  the  preassigned  error 

Qmnv»      The  search  for  point  P-,  is  called  the  E-Start 
max  r      1  

Problem, 
(ii)  Another  point  {x,y|  =  {xel,  yeel\  lies  on  tne  curve 
of  the  arctangent  function  such  that  x  =  x  -[_•   From 
here  point  {x,y}  starts  traveling  upwards  to  the  right. 
At  the  same  time,  the  position  of  the  line  L,  deter- 
mined by  the  points  {xel,  ^el^  and  {x'y}  »  is  GOm~ 
pared  with  the  tangent  line  at  \x,yj  of  the  arctan- 
gent function.   The  angle  between  these  two  lines  de- 
creases monotonically  with  x.   Point  {x,y|  stops  as 
soon  as  the  angle  $   approaches  zero,  where  line  L  co- 
incides with  the  full  tangent  line  tn.   The  point  of 
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contact,  t-j,  is  denoted  by  T^  =  (x^,,  y^i)  •  ^iG 
search  for  point  T-,  is  called  the  T-Problerfi  (Fig.  li^)  . 
(iii)  The  third  step  of  this  procedure  is  called  the  El- 
Problem.   While  the  E-Start  Problem  began  at  point 
(0,0{,  the  E-Problem  starts  at  the  point  of  tangency, 
T]_  =  (x-tl>  ytll  *    calculated  now.   The  point  \x,y) 
stops  at  point  P2  =  {*e2>    Je2 )    '      If  the  Procedure 
is  terminated  here,  it  is  called  a  1-Trip,  because  a 
polygonal  approximation  uses  one  full  tangent  line. 
The  horizontal  line  L  with  the  equation  y  =  ye2  is 
parallel  to  the  asymptotic  line  t^and  approximates 
arctan(x)  for  x  >  xq2. 
If  a  polygonal  approximation  with  two  full  tangent  lines 
is  desired,  the  procedure  of  the  T-Problem  (ii)  for  (xe2,  ye2] 
has  to  be  iterated  in  order  to  get  to  the  tangent  point 
^x,2,  yt2l  •   Prom  here,  the  procedure  of  the  E-Problerc  (iii) 
yields  point  ^x  o,  yeo}  .   This  polygonal  approximation  is 
called  a  2-Trip.   For  arguments  x  >  xe^,  the  line  L  with  the 
equation  y  =  yeo  is  used. 

If  E-Start  Problem,  T-Problem,  and  E-Problem  are  executed 
in  the  sequence 

E-Start,  T,  E,  T,  .  .  . ,  T,  E,  .  .  . 
where  E-  and  T-Problem  occur  N-times,  the  resulting  polygonal 
approximation  is  called  an  N-Trip. 

The  same  train  of  ideas  holds  for  all  N-Trips.   Hence  the 
"Upper  Chebyshev  Polygonal  Approximation  with  Zero  Boundary 
Error"  could  be  illustrated  by  the  2-Trip.   To  describe  N-Trips, 
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T={xt>ytl 


P={xe,yej 


Fig. 14.     The    T-  Problem. 
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the  following  symbols  are  used  for: 
Tangent  lines 

t0>    tl>    t2'    •*•>  %>  too    > 
corresponding  points  of  tangency 

m     rn     m  mm. 

Lq,    i]_,  i-2'    "'>    lY\>    lo°     ' 

x-coordinates   of   points   Tq,    -->,    m<>o 

xt0'  xtl>  xt2>  •"'  xtN'  xt  ; 
y-coordinates  of  points  Tq,    .  .  .  ,  T^ 

^t0>  ^tl'  yt2>  ••"  ^tN>  H    ; 
end  points  of  approximating  tangent  line  pieces  are 

P0,  ?lf    P2,  ...,  PN,  PN+1,  Poo  ,  where  PQ   =  TQ 

and  Poo  =  Tpo  as  tangent  lines  tQ  and  t00   are 

half  tangents; 

x-coordinates  of  points  P-j,  ••♦,  ?©<> 

xel>  xe2>  xe3>     '"'    xeN'  xe(N+l)  ; 
y-coordinates  of  points  P-p  ••-,    Poo 

yel>  Je2>    ye3>  •"'  yeN'  ye(N+l)  ' 
approximation  intervals  on  x-axis 

Iq*  I]_*  *2>  '  *  '  *    %'  ^-oo  ' 

Symbols  with  the  same  subscript  are  related  together,  e.g., 
tangent  line  to  contacts  the  graph  of  the  arctangent  function  at 
T,  which  has  the  coordinates  x, ^  and  y^o-   The  approximating 
line  piece  of  to  has  the  end  points  FS  and  Pl  whose  coordinates 

are  (xe3'  7e3)   and  {xel+'  Jek\    '      xe3  and  xel±   are  the  end  Points 
of  the  interval  It.   As  different  N-Trips  are  independent  from 

each  other,  there  is  no  need  for  left  subscripts.   N  has  a 

similar  function  as  the  k  of  the  k""1  approximation  stage. 


23 


The  2-Trip  mentioned  above  uses  the  tangent  line  tQ  at 
x  =  0,  but  not  necessarily  the  asymptotic  line  y  =  te/2  as  ap- 
proximating tangent  line.   The  eventual  result  of  the  2-Trip 
for  emax  =  0.01  is  yQ^  =  0.819  ....   As  ye^  <k/2,    the  value 
of  emax  has  to  be  increased.   If  emax  =  0.02,  ye2  =  1.015  ... 
emax  is  increased  until  ye^  =  ft/2.   This  particular  problem  is 
called  Pi-half  problem.   The  ensuing  approximation  is  called 
Asymptotic  2-Trip.   In  order  to  obtain  the  Asymptotic  N-Trip, 
the  same  reasoning  holds  for  any  other  N-Trip,  provided  N  is  a 
nonnegative  integer. 

Upper,  Lower,  and  Intermediate  Approximation 

As  the  arctangent  function  is  a  concave  function  in  the 
interval  (0,  ©o ) ,  the  upper  polygonal  approximation  e  ( x)  is  a 
tangent  approximation.   By  the  same  reason,  a  lower  approxima- 
tion g^(x)  of  arctan  ( x)  is  a  secant  approximation  in  (0,  oo   )• 
If  a  lower  and  an  upper  polygonal  approximation  have  the  same 
approximation  intervals  (Fig.  l£)  and  if  the  maximum  errors  of 
g. ( x)  and  gu(x)  are  uniform,  their  maximum  errors  are  equal. 
The  intermediate  approximation  of  a  lower  and  an  upper  approxi- 
mation with  common  approximation  intervals  is  defined  as 

1 
gj_(x)  =  -  (gu(x)  +  gj(x))  . 

1 

The  maximum  error  of  g^(x)  is  —  .  ©max«   Therefore  the  intermedi- 
ate approximation  allows  to  reduce  the  number  of  tangent  lines 
for  a  uniform  polygonal  approximation  with  a  preassigned 


2^ 


Fig.  15  .    Upper    and    lower    approximations  gu(x)  and  g,(x) 
of  function  f(x). 
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Fig.  16  .     Intermediate    polygonal     approximation    g,(x) 
of  function   'f(x) . 
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accuracy  (Fig.  16) . 

The  Step-division  Process 

The  E-,  E-Start,  and  T-Problems  still  to  be  solved  in  the 
next  subchapter  have  a  common  feature.   They  can  be  reduced  to 
finding  the  zero  crossing  of  a  test  function  h(x).   This  can  be 
done  by  establishing  an  appropriate  algorithm.   Subsequently, 
the  Step-division  Process  is  treated  in  detail. 

Let  h(x)  be  a  monotonic  and  continuous  function  defined  in 
a  closed  interval  (a,b).   It  is  assumed  that  there  are  two  argu- 
ments, xa  and  xb  in  (a,b)  whose  values  h(xQ)  and  h(xb)  have 
opposite  signs.   These  assumptions  necessarily  entail  a  zero 
crossing  for  h(x)  in  the  interval  (x  , x^).   The  Step-division 
Process  (SDP)  whose  purpose  is  to  locate  this  zero  crossing  can 
be  started  at  xQ .   There  are  two  types  of  Step-division  Processes. 

In  the  Unilateral  SDP,  x  approaches  the  zero  crossing  always 
from  one  direction.   The  process  starts  at  Xq  =  xa,  then  Xq  in- 
creases by  aQ  such  that  n  •  a$   =   x-^  -  xa,  where  n  is  a  positive 
integer.   After  computation  of  h(x-j_),  where  x-]_  =  xa  +  &q,    there 
are  two  cases  to  distinguish. 

Case  (a).   sign  (h(x1))  =  sign  (h(xQ)) 

This  means  that  there  is  no  zero  crossing  in  the  interval 
(xq,x-j)  and  that  the  argument  increases  again  such  that  X2  = 
X]_  +  a0.   This  process  is  continued  until  Case  (b)  occurs. 
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Case  (b) .   sign  (h(xi+1))  ^  sign  (h(xQ)) 
at  the  (i+1)    iteration.   Hence  the  zero  crossing  lies  in  the 
interval  (x.,  x.,-,).   In  order  to  bring  the  argument  back  into 
the  new  starting  position  xQ,  xi+i    is  decreased  by  aQ,  ( xQ  = 

X-+1  "  arJ •   Here  the  length  of  a^  is  divided  by  a  positive  in- 

a0 
teger  t,  (a-,  =  — )  such  that  t  >  1.   The  same  procedure,  de- 

t 
scribed  so  far,  is  repeated  for  a-^  and  the  new  value  of  Xq. 

Figure  17  shows  the  flow  chart  of  the  "Unilateral"  algo- 
rithm.  The  substitution  operation  that  is  denoted  by  a  reverted 
arrow  "^ —  "  avoids  subscripting  of  argument  x  and  step  a. 
a  < — b  is  read:   The  value  of  a  is  replaced  by  the  value  of  b. 
The  (a  -  5) -decision-box  terminates  computation  as  soon  as  the 
desired  accuracy,  indicated  by  6,  is  attained. 

In  the  Oscillatory  SDP,  x  approaches  the  zero  crossing 
alternately  from  either  direction.   As  in  the  "Unilateral  SDP", 
the  process  starts  at  x0  =  xa,  too,  and  x0  is  increased  like- 
wise.  The  value  h(x^)  is  calculated  for  each  argument  x^. 
Again,  there  are  two  cases  to  distinguish. 

Case  (a).   sign  (h(xi))  =  sign  (h(xi_1)) 
Case  (b).   sign  (h(xi))  ^  sign  (h(xi_1)). 
It  should  be  obvious  that  the  signs  of  two  succeeding  arguments 
xi-l  anc^  xi  are  comPare(3  here.   In  Case  (a),  x^  is  again  in- 
creased by  aQ.   In  Case  (b),  however,  there  exists  a  zero  cross- 
ing in  the  interval  (x^,  x^+1) .  x±+±   ~   xq   is  the  new  initial 

a0 
argument.   Step  (a)  is  divided  by  t,  (an  =  — ),  t  being  a  posi- 

t 

tive  integer,  greater  than  1.   The  new  step  a2  is  subtracted 

from  the  new  initial  argument  Xq,  and  the  same  process  is 
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USDP 


a  «•—  a/t 

. 

. 

x  ««—    x   -    a 

y 


Pig.  17.   Flow  chart  of  the  "Unilateral  step 
division  process"  (USDP) . 
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repeated  over  and  over,  resulting  in 

_  al  _  a0        a2  _  a0 
2    t     t2     3    t     t3 
Figure  18  shows  the  flow  chart  of  the  Oscillatory  SDP. 

The  Bisecting  Step  Process 

The  most  important  assignment  for  t  is  t  =  2.   In  this 
case,  the  interval  with  the  zero  crossing  is  bisected.   The 
features  of  both,  Unilateral  and  Oscillatory  SDP's,  can  be  com- 
bined.  Figure  19  shows  the  flow  chart  of  this  process,  called 
Bisecting  Step  Process.   The  advantage  of  this  SPD  over  pro- 
cesses with  t  >  2  is  due  to  the  fact  that  after  each  calculation 
of  y  =  f(x),  the  step  of  length  a  is  halved.   After  100  itera- 
tions, if  the  process  starts  with  an  initial  step  a  =  1.0,  the 

1 
error  of  the  approximation  is  at  most  .   For  the  processes 

2ioo 

with  t  >  2,  the  upper  bound  for  the  error  is 
t-100/t 

This  holds  for  the  worst  case  where  t  iterations  are  done  be- 
tween two  step  divisions.   For  the  specific  case  t  =  10,  the 
error  decreases  after  M  iterations  to 

1  1 


10: 


X    loO.LM    (1.259. ..)M 

where       100.1.M  =  100.1.M  =  (1.259. ..)M 

The  upper  bound  of  the  error  decreases  much  slower  than  for 
t  =  2: 
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^p   ^a 


IVY 


x,  •» —  xn  + 


ys  =   f(xs) 


i<*--  a/t 


i«-  +  a/t 


Pig.  18.   Flow  chart  of  the  oscillatory  step 
division  process  (OSDP). 


BSP 


A     , 

r     * 

x«* x  +    a 

1 

» 

y  =  h(x) 

Fig.  19.   Plow  chart  of  the  bisecting 
step  process  (BSP) . 
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1 
5max  =  —^ 

as  the  number  M  of  iterations  grows.  Although  in  case  t  =  10 

only  five  iterations  between  two  step  divisions  are  executed  on 
the  average 

1  1 


10:   6W 


max  10o.2.M   (i.585...)M 
the  SDP  for  t  =  10  is  still  slower  than  the  SDP  for  t  =  2. 

The  SDP  for  t  =  2  is  called  "BIT-BY-BIT  METHOD",  if  the 
numbers  used  in  the  process  are  represented  in  binary  notation. 
Likewise,  the  SDP  for  t  =  10  is  called  "DIT-BY-DIT  METHOD",  if 
the  numbers  used  in  the  process  are  represented  in  decimal  no- 
tation.  "DIT"  is  a  contraction  of  the  term  "decimal  digit" . 

In  the  "BIT-BY-BIT  METHOD"  as  well  as  in  the  "DIT-BY-DIT  METHOD", 

a  1 
the  division  —  =  —  at  each  step  change  takes  very  little  time  to 

t  t 
be  performed.   Suppose  "BIT-BY-BIT"  (BBB)  and  "DIT-BY-DIT"  (DBD) 

are  started  with  unit  step  length,  the  following  sequences  of 

steps  ooour. 

BBB:   a  =  1.(2)*  °'1{2)>    °-01(2)'  0.001(2)>  ... 

DBD:   a  =  1.    ,0.1    ,0.01    ,  0.001    ,  ... 
The  subscript  "(2)"  indicates  binary  notation.   For  all  step 
lengths,  the  mantissa  is  one  digit  long,  whereas  in  other  SDP '  s 
the  length  of  the  mantissa  of  the  number  representing  step 
length  is  increased  by  each  division.   Thus  computation  time  is 
larger,  and  the  mantissas  of  the  step  lengths  are  truncated. 
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The  Newton  Process 

In  the  "Newton  Process"  (NP),  the  conditions  to  be  satis- 
fied by  h(x)  are  the  same  as  those  for  the  "Step-division  Pro- 
cess", with  the  exception  that  there  exists  for  h(x)  a  continuous 
monotonic  first  derivative  h'(x)  in  the  interval  (x  ,  x,) . 
These  assumptions  are  sufficient  for  the  convergence  of  the  NP. 

There  are  functions  h(x)  with  non-monotonic  continuous  first 
derivative.   h(x)  with  an  inflection  point  at  the  zero  crossing 
shows  special  interest.   There  are  two  cases  to  distinguish. 

(a)  If  h(x)  is  concave  for  arguments  less  than  the  argu- 
ment of  zero  crossing  and  convex  for  arguments 
greater  than  the  argument  of  zero  crossing,  the  NP 
is  convergent. 

(b)  If  h(x)  is  convex  for  arguments  less  than  the  argu- 
ment of  zero  crossing  and  concave  for  arguments 
greater  than  the  argument  of  zero  crossing,  the  NP 
is  either  convergent  or  divergent,  depending  on  the 
local  properties  of  h(x). 

The  assertion  of  case  (a)  is  thoroughly  treated  in  Appendix  B. 

The  flow  chart  of  the  "Newton  Process"  is  shown  in  Fig,  20. 

The  E-Start  Algorithm 

The  N-TRIP  problem  starts  with  finding  the  argument  xel 
where  the  maximal  error  emax  of  the  half  tangent  line  tQ  occurs. 
The  input  data  for  this  algorithm  are: 
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NP 


V 

f   =   f(x) 
f«    -  ft  (x) 

Fig.  20.   Flow  chart  of  the  Newton  process  (NP) 
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xto  =  °>   yto  =  0;   yto  =  x  • 

The  "E-Start  Algorithm"  (E-START)  searches  for  the  zero  crossing 
of  the  test  function 

h(x)  =  (x  -  arctan(x))  -  emax  . 
As  initial  step,  a  =  1.0  is  appropriate.   Figure  21  shows  the 
flow  chart  of  the  E-Start  Algorithm  that  uses  a  DBD  method.   Ex- 
pressions in  the  different  blocks  are  written  in  FORTRAN  lan- 
guage.  The  "equals"  signs  in  the  rectangular  assertion  boxes 
refer  to  "substitution",  that  were  denoted  by  "  -« —  "  in  non- 
FORTRAN  flow  charts.   Table  1  explains  symbols  used  in  FORTRAN 
language.   The  output  data  of  the  E-START  are  the  approximate 
values  for  xQl  and  ye}_. 

Experience  with  the  IBM  1620  electronic  digital  computer 
has  shown  that  the  E-START  computation  time  can  be  reduced  by  a 
factor  3  when  using  a  Newton  Process.   The  test  function 

h(x)  =  (x  -  arctan(x))  -  emax 
satisfies  all  the  requirements  for  a  NP.   The  calculation 
1 


h'(x)  -  1 


1  +  x2 

leads  to  the  recurrence  relation 
h(x) 


xn+l  ~  *n  " 


xn+l 


h'(x) 

( x  -  arctan( x) )  -  e„ 


1  +  x2 
Figure  22  shows  the  FORTRAN  flow  chart  of  the  E-START  that  uses 
an  NP. 
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E-START 


XE  =  A 


XE  =  XE+A 


H> 


|YEE  =  ATAN(XE)  | 

1 EPSXE  =  (XE-YEE)-EMAX 


h    >0 


.  l-::-A 


( ABS (EPSXE) -DELTA )- 


XE=XE+AI 


XE  =  XE-AI 


J- 


YEE  =  ATAN(XE) 


_L 


|  EPNXE  =  (XE-YEE)-EMAX| 


EPNXE-::-EPSXEJ 


AI 


.  l-::-AI 


>0 


EPSXE  =  EPNXE 


YE  =  XE 


1  YE  I  =  .£-»•(  YE +YEE)| 


Fig.  21.   FORTRAN  flow  chart  of  the  E-START 
algorithm  with  DBD  method. 
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Table  1.   Explanation  of  FORTRAN  symbols. 


FORTRAN 

variable 

symbol 


Standard 
notation 


Name  of  variable  that 
approaches  . . . 


XE 

YE 

YEE 

YEI 

EPSXE 

EPNXE 

EPS 


y->ye 
y->yel 


3(X) 


.  argument  xQ 

.  y 
J  ee 

*  ^ei 
.  DELTA 

.  DELTA 

rror 


XT 

YT 

YDT 

B 

EPSXT 

EPNXT 


y-*yt 
y— *yt 
b-*bt 

h(xp) 
h(xa) 


...  x+_  )  coordinates  of 

t  ) 
...  y+  )  point  of  tangency  Tn 

...  y-t    slope  of  tangent  line 

. . .  bt   y-intercept  of  "   " 

error  at  argument  x 

error  at  argument  x  +  h 


A 

Al 

DELTA 

EMAX 
ATAN 
ABS 


a0 

a1,  a2, 


max 
arctan 


constant  step 

variable  step 

constant  reference  error  that 

determines  arithmetic  accuracy 

maximal  error 

arctangent  function 

absolute  value 
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E-START 


r11 

|XE   = 


YEE   =   ATAN(XE) 


J_ 


IePSXE    -    (XE-YEE)-EMAX 

i  

|  XE   =   XE+a[«<0  (ePSXe) 


1 


>0 


|YDE   =   l./(l.+XE-::--::-2.) 


i 


(aBS  ( EPSXE  )-DELTa) 

>°  

|XE    =  XE-EPSXE/(1.0-YDE) 

I 

| YEE    =   ATAN(XE)    | 
[EPSXE    =    (XE-YEE)-EMAX 

1 


I  YE    =   XE 

I 


|  YE  I    =    .5---(YE+YEE)  | 


Fig.    22.      FORTRAN    flow   chart    of    the   E-START 
algorithm  with  Newton   process. 


38 


The  T-Algorithm 

In  the  subchapter  on  "Upper  Chebyshev  Polygonal  Approxima- 
tion with  Zero  Boundary  Error"  it  was  explained  how,  at  the  point 
of  tangency,  the  angle  0   becomes  zero.   There,  Pig.  llj.  makes 
plausible  that  the  length  of  the  directed  segment   FP   tends 
toward  zero,  as  x  tends  toward  x^ .   This  yields  the  test  function 

h(x)  =  PPn  =  (y  -  ye)  -  y(x  -  xe) 
or 

h(x)  =  (arctan(x)  -  jQ)    -  ( ( x-xQ) /( 1+x2);  . 

Input  data  are  xQ  and  ye,  obtained  from  a  previously  computed 
"E-START  Algorithm"  or  an  "E-Algorithm" .   Output  data  are  xt  and 
yt.   Figure  23  shows  the  "T-Algorithm"  solved  with  the  DBD 
Method. 

The  function  h(x)  meets  requirements  for  the  "Newton  Pro- 
cess".  The  first  derivative  of  h(x)  is: 

(x  -  xQ)2x 
h'(x)  =  — —  =    (x  -  xe)2xy^   . 

(1  +  x^ ) d 

Functions  h(x)  and  h'(x)  yield  the  following  recurrence  relation" 

(y  -  ye)  -  J(x  -  xQ) 
xn+l  ~  xn  " 

where  xQ  =  xQ .   Figure  2\\.   shows  the  FORTRAN  flow  chart  of  the 
T-Algorithm  with  Newton  Process. 
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XT    =   XE+A 


YDT    =    l./(l-+XT-::--::-2.) 

I 

lYT   =   ATAN(XT)     I 

*  ^ 

IEPSXT   =    (YT-YE)-YDT-::-(XT-XE) 


<0  , 


|  XT    =  XT+AjO — (  EPSXT) 

j  >o 


Al   =    .1-x-A 


_L 


|ABS( EPSXT) -DELTA 

^  >0 
<Or{EPSXT  )-[    >0 


XE   =  XE+A1 


XE   =  XE-A1 


1,11 

|YDT  =  l./(l.+XT-*-*2.)  | 

# 

|YT  =  ATAN(XT) 

♦ 

|  EPNXT  =  (YT-YE)-YDT-::-(XT-XE)  | 

* 

J 

(EPNXT-::-EPSXT  ) 
fr^O    |  >0 

^  1 

| EPSXT  =  EPNXT 

_i 

=  YT-YDT-"-XT   r*" 


UT 


Pig.    23.      FORTRAN   flow   chart    of   the    T-Problem 
with  DBD  method. 
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I  XT  =  XE+A 


| XT  =  XT+A  [— -(  EPS 


|YDT  =  l./(l.+XT -::•-::- 2.  )| 

|YT  =  ATAN(XT)  | 

I  

JEPSXT  =  ( YT -YE  )-YDT-*(  XT-XE) 
\  n 


~1     , 

(  ABS(EPSXT) -DELTA) 

I    ° 

[XT  =  XT-EPSXT/(  2  .  -::-XT-::- ( XT-XE )  -::- ( YDT-*-«-2  .  )  )  I 

1 

lYDT   =    l./(l.+XT-*-*2.)    I 

I  , 


YT   =   ATAN(XT) 


EPSXT    =    (YT-YE) -YDT-"-( XT-XE)    | 


J 


£ 


B   =  YT-YDT-::-XT 


IT 


Fig.  21;.   FORTRAN  flow  chart  of  the  T-Algorithm 
with  Newton  process. 
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The  E-Algorithm 

For  an  argument  x,  the  error  between  the  approximating 
tangent  line  and  arctan(x)  is: 

e(x)  =  (mx  +  b)  -  arctan(x). 
The  aim  of  the  E-problem  is  to  find  the  argument  x  =  xe  for 
which  e(xe)  =  emax  holds,  and  therefore  the  E-problem  reduces 
to  finding  a  zero  of  the  test  function 

h(x)  =  ((mx  +  b)  -  arctan(x))  -  emax  . 
Input  data  for  the  "E-Algorithm"  are  x^  and  j^,    obtained  from 
the  previously  computed  T-Algorithm.   Output  data  of  the 
"T-Algorithm"  are  xQ  and  yQ .   Figure  25  shows  the  FORTRAN  flow 
chart  of  the  "E-Algorithm"  with  DBD  Method. 

Likewise,  h(x)  meets  requirements  for  use  of  the  "Newton 
Process".   The  first  derivative  of  h(x)  is 

h'(x)  =  m  -  (1/(1  +x2))  . 
The  recurrence  relation  reads: 


(m^  +  b)  -  arctan(xn) 
m  -   1/(1  +  xn^) 


E-Algorithm  with  Newton  Process. 

The  N-TRIP 

The  E-START,  the  E-  and  the  T-algorithms  are  the  building 
blocks  for  an  algorithm  to  solve  the  N-TRIP  problem.   The  5-TRIP 
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XE   =  XT+A 


YE    =  YDT-::-XE+B 


_ 


YEE    =   ATAN(XE) 


X 


EPS   =  YE -YEE 


_L 


EPSXE    =  EPS-EMAX 


|  XE   =  XE+A  \<  *°(ePSXE  ) 

ir  >o 


Al   =    .1*-A 


(ABS( EPSXE)  -DELTA  )- 
<0H  EPSXE      H   >0 
XE    =   XE+All     |XE    =   XE-A1 


*=0 


_ 


YE  =  YDT-::-XE+B 

H2 1 

YEE  =  ATAN(XE) 

1 1 l 

[  EPS  =  YE -YEE   | 

I  EPNXE  =  EPS-EMAX 


[epnxe-»-epsxe  ) 


^<Q 


Al 


■  1-*A1| 


_ 


>0 


JEPSXE  =  EPNXE 


Pig.  25.   FORTRAN  flow  chart  of  the  E-Algorithm 
with  DBD  Method. 
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XE  =  XT+A 


XE  =  XE+A 


£0 


9   ? 

YEE  =  ATAN(XE)1 

j         . 

EPSXE  =  ( (YDT*-XE+B)-YEE)-EMAX  | 

epsxeJ 


1 


>0 


saO 


|YDE  =  l./(l.+XE-::--::-2.) 

(ABS( EPSXE) -DELTA  ) 

|XE  =  XE-EPSXE/(YDT-YDE)  I 

I 

IYEE  =  ATAN(XE)  1 

I  EPSXE  =  (  (YDT-::-XE+B)-YEE)-EMAX  | 


* 

|ye  = 

YDT-: 

.-XE+B  | 

1 

|yei 

=  .9- 

>-(YE+YEE)   1 

Fig.  26.   FORTRAN  flow  chart  of  the  E-Algorithm 
with  Newton  process. 
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with  emax  =0.01  serves  as  an  example.   First,  the  E-START  prob- 
lem is  solved,  yielding  {xei>yei'l>  followed  by  the  T  problem, 
yielding  jxtl^tl}*   Recurrent  solving  of  the  last  two  problems 
leads  to  the  sequence: 

{xt2,yt2},  (xe2,ye2},  {xt3,yt3},  {*tk'Ytk\' 

lxeVye^'  (Xt5'yt5l'  iXe5'ye5i- 

Figure  27  shows  the  required  FORTRAN  flow  chart. 

The  Asymptotic  N-TRIP  Algorithm 

When  yeN+1  =  V2>  tne  N-TRIP  becomes  an  Asymptotic  N-TRIP. 
For  the  N-TRIP,  yeN  is  a  function  of  emax.  Henceforth  the  test 
function  reads: 

h(x)  =  yeN+1(emax)  =  U/2) 
This  test  function  is  transcendental  and  satisfies  the  require- 
ments for  the  Newton  Process.   Being  a  composite  transcendental 
function,  it  is  awkward  to  transform  into  a  closed  expression 
that  would  be  needed  for  the  speedy  Newton  Process.   In  this  case 
the  DBD  Method  proves  to  be  more  workable.   It  starts  with  a  suf- 
ficiently small  value  for  e    for  which  h(x)  is  negative.   In 
case  of  a  2-TRIP,  emax  =  0.01  is  sufficiently  small.   emax  grows 
gradually  by  increments  of  0.01  until  h(x)  becomes  negative. 
The  procedure  is  then  continued  in  the  usual  way  (Figs.  28  and 
29) .   The  algorithm  to  find  the  zero  of  the  above  test  function 
is  called  pi-half  algorithm. 
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Pig.    27.      FORTRAN   flow   chart    of  an  N-TRIP. 
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NTRIP 
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DEL 


EMAX 

=  EMAX-DELT 

DELT 

=  0.1-::- DELT 

[emax 

=  EMAX+DELT 

EMAX  =  EMAX+DELT 


Fig.  28.   FORTRAN  flow  chart  of  the  asymptotic 
N-TRIP  with  unilateral  SDP. 
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* 
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■    .  l-::-DELT  | 
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I 

EMAX 
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' 
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I 

■ 

F 

Fig.  29.   Flow  chart  of  asymptotic  N-TRIP 
with  oscillatory  SDP. 
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Tables  2  through  11  show  the  results  of  the  Asymptotic  1-, 

2-,  .  .  .,  10-TRIP.   The  way  to  use  these  tables  is  indicated 

for  the  Asymptotic  1|-TRIP  of  Table  12.   emax  =  .1^01706E-01  is 

the  value  of  the  maximal  error  of  the  Asymptotic  lj_-TRIP,  and 

"E-01"  means  10"   in  standard  mathematical  notation.   The  first 

quadruple  of  the  numbers  XE,  YE,  YEE,  YEI,  called  E-START 

values,  reads  in  standard  notation  x  ,,  ^ e\>    ^eel'  yeil'   Tiie7 

are  the  result  (output  data)  of  the  E-START  algorithm  and  are 

used  as  input  data  for  the  ensuing  T  algorithm,  whose  output 

data  are  presented  in  the  next  quadruple  of  the  numbers  XT,  YT, 

YDT,  B.   In  standard  notation,  they  are  indicated  by  xtl,  y^l' 

y  ,  ,  m,  .,  .   The  next  quadruple  of  E-values  has  subscript  2;  xQ2> 

J   of    J      o>    Y    -o'      T*16  succeeding  quadruple  of  E-values  carries 
e  c.        e  e  c.        e  i  <i 

again  subscript  2.   With  each  new  quadruple  of  E-values,  the 
subscript  increases  by  one.   The  last  quadruple  of  numbers  is 

e5'   ©5*   ee5'  ^ei5* 

Tables  2  through  12  show  the  values  of  emax  for  different 

Asymptotic  N-TRIPS.   M  denotes  the  total  number  of  tangent  lines 
of  the  polygonal  approximation  for  the  interval  (0,  oo  )  .   N  is 
the  number  of  full  tangents  in  (0,  oo)  (Table  13) . 

In  the  chapter  on  Quasi-equal  Error  Method  it  was  recog- 
nized that  the  upper  polygonal  approximation  for  the  second 
stage  has  maximum  errors  with  uniform  Cheybshev  norm.   This  ap- 
proximation is  identical  with  the  upper  polygonal  approximation 
of  the  Asymptotic  1-Trip.   This  assertion  is  inferred  from  the 
observation  that  both  approximations  employ  the  half  tangent 
lines  at  x  =  0  and  x  =  00  and  use  the  same  number  of  full  tan- 
gents.  If  the  uniqueness  of  the  Asymptotic  1-Trip  approximation 
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is  taken  for  granted,  this  fact  implies  the  above  identity. 

The  fact  that  the  upper  polygonal  approximation  of  the 
second  stage  of  the  Quasi-equal  Error  Method  has  uniform  Cheby- 
shev  norm,  suggests  to  ask  if  the  upper  polygonal  approximation 
of  the  third  stage  has  maximum  errors  with  uniform  Chebyshev 
norm,  too.   This  approximation  corresponds  to  the  Asymptotic 
3-Trip.   Similar  reasoning  as  above  leads  to  the  assertion, 
"If  the  upper  polygonal  approximation  of  the  third  stage  of  the 
Quasi-equal  Error  Method  has  maximum  error  with  uniform  Cheby- 
shev norm,  then  it  is  identical  with  the  upper  polygonal  ap- 
proximation of  the  Asymptotic  3-Trip."   Since  the  full  tangent 
line  pt-i  of  the  first  stage  approximation  of  the  Quasi-equal 
Error  Method  is  used  as  approximant  in  all  other  higher  stages, 
the  argument  2xtl  °^  ^s  P°int  °f  tangency  is  identical  with 
oX,p  (Pig.  9) .   This  observation  together  with  the  two  assertions 
mentioned  above  lead  to  the  conclusion  that  argument  x,-,  of  the 
full  tangent  line  belonging  to  the  Asymptotic  1-Trip,  is  equal 
to  the  argument  xtp  °f  the  middle  tangent  of  the  three  full  tan- 
gent lines  of  the  Asymptotic  3-Trip.   Approximate  values  for 
these  arguments  can  be  found  in  Tables  2  and  l±,    where 

xtl  =  2.3878... 
and  xt2  =  2.6311. . . 

These  two  real  numbers  are  different.   Therefore  it  is  proven 
that  the  upper  polygonal  approximation  of  the  third  stage  of  the 
Quasi-equal  Error  Method  has  maximum  error  with  uniform  Chevy- 
shev  norm  must  be  negated. 

Similar  reasoning  for  the  upper  polygonal  approximation  at 
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the  fourth  stage  of  the  Quasi-equal  Error  Method  leads  to  the 
conclusion  that  this  approximation  corresponds  to  the  Asymptotic 
7-Trip  and  that  x. -,  must  be  equal  to  argument  x  ,  of  the  Asymp- 
totic 7-Trip.   Table  8  supplies  the  value  of  xt>  =  2. 82^63.   As 
xtl  -=f   x^l,  it  is  concluded  that  the  above  mentioned  conjecture 
has  to  be  negated  for  the  fourth  stage  polygonal  approximation 
of  the  Quasi-equal  Error  Method,  too.   A  general  negation  proof 
is  not  given  here.   I  conjecture  such  a  proof  would  show  that 
the  progression  xtl(l_Trip),  xt2(3_Trip)>  ^ 7-Trip) >  ■  '  " 
is  monotonic. 


Table   2.      Asymptotic    1-TRIP. 
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MAX  = 

.195&52E-00 

XE   = 
YEE   = 

.96131J.3E  00 
.7&5691E   00 

XT   = 
YDT    = 

.238783E  01 

.1IJ.9215E  00 

XE   = 
YEE   = 

.50I4.573E  01 
.137511^  01 

XEN   = 

•  50IJ.573E   01 

YE   = 
YEI    = 

.96l3i^3E 
. 863517E 

00 

00 

YT    = 
B   = 

.11714-20E 
. 817596E 

01 

00 

YE   = 
YEI   = 

.157080E 
.1J+7297E 

01 
01 

YEN   = 

.157080E 

01 

Table   3.      Asymptotic    2-TRIP. 


MAX   = 

.991702E-01 

XE    = 
YEE    = 

.72927I4-E  00 

.63010l|E   00 

XT   = 
YDT   = 

.  1i4.2J4.75E  01 
.330JOI4.3E  00 

XE    = 
YEE   = 

.2369I4.7E  01 
.1171W  01 

XT   = 
YDT    = 

.J4.9586 IE  01 
.390811E-01 

XE    = 
YEE   = 

.100506E   02 

.H1.7163E  01 

XEN   = 

.100506E   02 

YE 
YEI 

- 

.72927I4E  00 
.679689E  00 

YT 

B 

= 

.958811E  00 
.^88582E  00 

YE 
YEI 

= 

.127061E  01 

.122102E   01 

YT 
B 

= 

.137180E  01 
.117801E  01 

YE 
YEI 

= 

.157080E  01 
.152121E  01 

YEN 

= 

.157080E  01 

Table   I4..      Asymptotic   3-TRIP. 
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MAX  = 

•599589E-01 

XE  = 
YEE  = 

.6011^27E  00 
.5Ij.llJ.68E  00 

XT  = 
YDT  = 

.107313E  01 
.J+64766E  00 

XE  = 
YEE  = 

.163520E  01 

.102193E  01 

XT  = 
YDT  = 

.263113E  01 
.126217E  00 

XE  = 
YEE  = 

.  1^.07 899E  01 

.133038E  01 

XT  = 
YDT  = 

.828891^  01 

.li|3i|-59E-01 

XE  = 

YEE  = 

.166581E  02 

.15108I4B  01 

XEN  = 

.166581E  02 

YE  = 
YEI  - 

.60j:Jj?7E  00 
.  "  Lj4  IE  00 

YT  = 

B  = 

.82066IE  00 

.32190I4E  00 

YE  = 
YEI  = 

.108189E  01 

.105191E  01 

YT  = 

B  = 

.120759E  01 
.875i98E  00 

YE  = 

YEI  = 

.13903I4E  01 
.136036E  01 

YT  = 

B  = 

.11^5073E  01 
.133182E  01 

YE  = 

YEI  = 

.157080E  01 

.1514-082E  01 

YEN  = 

.157080E  01 

Table  5.   Asymptotic  i^-TRIP. 
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MAX  = 

.1^0170oE-01 

XE   = 
YEE   = 

.5l8£1^0E  00 
.J4.78370E   00 

XT    = 
YDT    = 

.882930E  00 
•561935E   00 

XE   = 
YEE   = 

.128723E  01 
.910325E  00 

XT    = 
YDT    = 

.187333E   01 

.221760E  00 

XE    = 
YEE   = 

.262319E  01 

•120659E   01 

XT   = 

YDT   = 

,1+o£o6oe  01 

•  57ij468E-01 

XE    = 
YEE    = 

.616301E  01 

,li|.099i|E   01 

XT   = 
YDT    = 

,12l\13l\E  02 
.6I|4780E-02 

XE   = 

YEE    = 

.2i|.8803E  02 
.153062E   01 

XEN   = 

.21+8803E   02 

YE   = 

YEI    = 

.5i85i^0E  00 

.I^981|55E   00 

YT    = 
B   = 

.72330ijE    00 

.22715I4E  00 

YE   = 
YEI    = 

.950I4-96E  00 
.930I4.11E  00 

YT   = 
B   = 

.1080i|.7E  01 
.665039E   00 

YE   = 
YEI    = 

.12L1.676E  01 
.122667E  01 

YT   = 

B   = 

.132876E  01 

.109606E   01 

YE    = 
YEI    = 

.11+50  HE  01 
.ll|3002E    01 

YT   = 
B   = 

.lij.90ij.lE  01 
,liJ.1037E  01 

YE    = 
YEI    = 

.157080E  01 
.155071E  01 

YEN   = 

.157080E  01 

Table   6.      Asymptotic   5-TRIP. 


5k 


CMAX   = 

.287930E-01 

XE   = 
YEE   = 

.1+59609E    00 

.530816E  00 

XT    = 

YDT    = 

.760809E  00 
.633380E  00 

XE   = 
YEE   = 

.1080J±7E  01 
.  82I4.O6OE    00 

XT    = 
YDT    = 

.11+9278E  01 

.309751E    00 

XE   = 
YEE   = 

.I98I4.OOE  01 
.110393E   01 

XT    = 
YDT    = 

.275271E  01 
.116585E  00 

XE    = 
YEE    = 

.3757i4-8E  01 
.131069E   01 

XT   = 
YDT   = 

.571826E  01 
.2967I4.9E-OI 

XE    = 

YEE   = 

.8639ip.E  01 
.11^555^  01 

XT    = 
YDT   = 

.17314-13E  02 
.331U2E-02 

XE   = 
YEE   = 

.3V7210E  02 
.1514-200E    01 

XEN    = 

•3V7210E   02 

YE   = 
YEI   = 

.I4.59609E  00 

.IJ45212E   00 

YT   = 

B    = 

.650383E  00 

.168502E   00 

YE   = 
YEI    = 

.852853E  00 

.838i4.56E   00 

YT   = 
B   = 

.980566E  00 
.518175E  00 

YE    = 
YEI    = 

.113272E  01 
.111832E  01 

YT   = 
B   = 

.12223I|E  01 
.9011A7E   00 

YE    = 
YEI    = 

.1339i|8E   01 

.132509E  01 

YT    = 

B   = 

.139767E  01 

.122798E   01 

YE    = 

YEI    = 

.1i4.8J4.35E  01 
.1I4.6996E  01 

YT   = 
B   = 

.151319E  01 
.11,551 2E  01 

YE   - 
YEI    = 

.157080E  01 
.155614-0E   01 

YEN   = 

.157080E   01 

Table    7.      Asymptotic    6-TRIP. 
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UMAX  = 

. 21650 3E-01 

XE   = 
YEE   = 

.1J.15139E  00 
.  393J4-89E  00 

XT   = 
YDT    = 

.671+461^  00 

.687331E   00 

XE   = 
YEE   = 

.9I-I-1517E  00 
.755285E   00 

XT    = 
YDT   = 

.12606I4E  01 
.386217E  00 

XE   = 

YEE    = 

.162395E  01 
.101885E  01 

XT    = 
YDT   = 

.213058E  01 
.180526E  00 

XE    = 
YEE   = 

.27i|829E  01 
.122182E   01 

XT   = 
YDT    = 

.37i+352E  01 
.6660I+7E-OI 

XE    = 

YEE   = 

.505605E  01 
.137553E  01 

XT    = 
YDT    = 

.7%5W  01 
.168201E-01 

XE   = 

YEE   = 

.115114-6E  02 
.  li4.8ij.i7E    01 

XT   = 
YDT   = 

.230762E  02 
.187437E-02 

XE   = 
YEE   = 

.lj.6i8l3E  02 
.l5^91i4E   01 

XEN   = 

.ij.6l8l3E  02 

YE 

YEI 

= 

.I^l5l39E  00 
.ij.0J+311jE   00 

YT 
B 

= 

•593382E  00 
.129801E   00 

YE 
YEI 

= 

.776935E  00 
.766110E  00 

YT 
B 

= 

.900186E   00 

.I4.13305E  00 

YE 
YEI 

= 

.10i|.050E  01 
.102968E    01 

YT 
B 

= 

.113196E  01 
•  7i|7339E   00 

YE 
YEI 

= 

.12l^3l4-8E  01 
.123265E   01 

YT 
B 

= 

.130976E  01 

.  1060i|3E   01 

YE 
YEI 

I 

.139718E    01 

.138636E  01 

YT 
B 

= 

.  ii+J+o  71+E  01 

,13121ijE   01 

YE 
YEI 

= 

.150582E  01 
.lij.9ij.99E   01 

YT 
B 

= 

,l£27i|.9E  01 
.I^8j+23E   01 

YE 
YEI 

= 

.157080E  01 
.155997E   01 

YEN 

= 

.157080E   01 
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Asymptotic    7-TRIP. 


EMAX  = 

.168726E-01 

XE   = 
YEE   = 

.380155E  00 
.363283E  00 

XT    = 
YDT   = 

.609519E  00 
.729122E   00 

XE   = 
YEE   = 

.8I].06S9E  00 
.6990lj.6E   00 

XT   = 
YDT    = 

.1102l+i|E  01 
.!|.5l390E   00 

XE    = 

YEE   = 

.139126E  01 
.914-7582E   00 

XT   = 
YDT    = 

.1761jl|.2E  01 
.21|3122E    00 

XE    = 
YEE   = 

.219951E  01 
,lli|4o8E  01 

XT   = 
YDT   = 

.28256^  01 
.111306E  00 

XE    = 
YEE    = 

.359701E    01 

.129963E  01 

XT    = 
YDT   = 

.i^85707E  01 
.1+06650E-01 

XE   = 
YEE    = 

.652618E  01 

.lklQ7$E    01 

XT    = 
YDT   = 

.983695E  01 

.102285E-01 

XE   = 
YEE   = 

.11^79l6E   02 

.150329E  01 

XT   = 
YDT   = 

.296198E  02 
.113852E-02 

XE    = 
YEE   = 

.592620E  02 
•155392E   01 

XEN   = 

.592620E   02 

YE   = 
YEI    = 

.380155s  00 

.371719E    00 

YT    = 
B   = 

.5V7389E   00 

.102976E  00 

YE    = 
YEI    = 

.715919E  00 

.707I4-83E  00 

YT   = 
B   = 

.83I1085E  00 
.336I+54E  00 

YE    = 
YEI    = 

.96I1455E  00 
.956018E  00 

YT   = 
B   = 

.105518E  01 
.626209E  00 

YE   = 
YEI    = 

.116096E  01 
.115252E  01 

YT   = 

B   = 

.123065E  01 
.916138E  00 

YE   = 
YEI    = 

.1316S1E  01 
.130807E  01 

YT   = 

B   = 

.136775E  01 
.117023E  01 

YE   = 
YEI    = 

.ll;3562E  01 
.li|2719E   01 

YT   - 
B   = 

.llj.69ii.9E  01 
.136807E  01 

YE    = 
YEI    = 

.152016E  01 

.15H73E   01 

YT   = 
B   = 

.153705E   01 

.150332E  01 

YE   = 
YEI    = 

.157080E   01 

.156236E  01 

YEN   = 

.157080E  01 

Table  9.   Asymptotic  8-TRIP. 
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MAX   = 

.135193E-01 

XE    = 

YEE    = 

.351771E    00 

.338252E  00 

XT   = 
YDT   = 

•558519E   00 

.762228E  00 

XE   = 

YEE    = 

.763517E  00 
.652096E  00 

XT    = 
YDT   = 

.986675E  00 
.506707E  00 

XE   = 
YEE   = 

.122730E  01 
.887097E  00 

XT    = 
YDT    = 

.152177E  01 
.301587E  00 

XE   = 

YEE   = 

.18514.59E   01 

.107628E  01 

XT    = 
YDT    = 

.229891E  01 

.159109E   00 

XE    = 

.282282E    01 

YEE   = 

.123033E  01 

XT   = 
YDT   = 

.358872E  01 
.720517E-01 

XE    = 
YEE    = 

.1;53739E  01 
.135387E   01 

XT   = 
YDT   = 

.609852E  01 
.261835E-01 

XE   = 
YEE   = 

.817131^  01 

,li|4902E   01 

XT    = 

YDT   = 

.122951E   02 

.657156E-02 

YE   = 

.351771E    00 

YEI    = 

.3^0 12E   00 

YT    = 
B   = 

.509360E  00 

.836J4.1i|E-01 

YE   = 
YEI    = 

.665615E  00 
.658856E  00 

YT    = 
B   = 

.778691E  00 
.278736E  00 

YE   = 

YEI    = 

.900616E  00 

.893857E   00 

YT    = 
B   = 

.9891]-26E   00 

.  530I4.79E  00 

YE    = 
YEI    = 

.108980E  01 
.10830ljE   01 

YT   = 
B   = 

.  H60I4.9E  01 
.  79ij-717E  00 

YE    = 
YEI    = 

.12I4.385E  01 

.123709E   01 

YT    = 
B   = 

.12990i|E  01 
.lOJ4.Olj.6E   01 

YE   = 
YEI    = 

.136739E   01 

.136063E  01 

YT    = 
B   ■ 

.H4.0827E  01 

.12I|.859E   01 

YE   = 
YEI   = 

.lk.62$l±E   01 

.1I4.5578E  01 

YT    sb 
B   = 

.H4.896I4E  01 
.II4.O88I4E   01 

Table   9    (Concl.) . 
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XE    = 
YEE   = 

.18V717E   02 

.151671E  01 

XT    = 
YDT    = 

.369729E  02 
.730998E-03 

XE    = 
YEE   = 

.739638E  02 
.155728E   01 

XEN   = 

.739638E   02 

YE    = 
YEI    = 

.153023E 
.15231]-7E 

01 
01 

YT   = 

B   = 

.15^376E 

.151673E 

01 
01 

YE    = 
YEI    = 

.157080E 

.  1561|0I|E 

01 

01 

YEN  = 

.157080E 

01 

Table  10.   Asymptotic  9-TRIP. 
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EXAM 


XE 
YEE 


XT 
YDT 


XE 

YEE 


XT 
YDT 


XE 
YEE 


XT 
YDT 


XE 
YEE 


XT 
YDT 


XE 
YEE 


XT 
YDT 


XE 
YEE 


XT 
YDT 


XE 
YEE 


XT 
YDT 


.110755E-01 
.328190E  00 

.31711I4E  00 

.5l7l8i^E   00 
~~967E   00 


.7022I|.0E   00 
.612227E   00 

.897667E  00 
.553769E  00 

.HOl^E  01 

.835126E  00 

.  13i^8i3E  01 
.35i|932E   00 

.161696E  01 

.101692E  01 

.195710E  01 
. 207029E  00 

.23k607E  01 

.116787E  01 

.287^21+E  01 
.107977E  00 

.350081iE  01 

.129256E  01 

.kh.2k.9hE   01 
.lj.o,5906E-01 

.557310E  01 
.139325E  01 

.7V70£J±E  01 
.176026E-OI 


YE   = 
YEI    = 

.328190E  00 
.322652E  00 

YT    = 

B   = 

.k77300E  00 
.692587E-01 

YE   = 
YEI    = 

.623303E  00 
.617765E  00 

YT   = 
B   = 

.731525E  00 
,23i|I|.2i|E   00 

YE   = 
YEI    = 

.81|6202E  00 
.8IJ.O66I4E   00 

YT   = 
B   = 

.932583E  00 

.1+51+090E   00 

YE   = 
YEI    = 

.102800E    01 

.1022I4.6E  01 

YT   = 
B   = 

.10981+2E  01 
.6932i+2E   00 

YE   = 
YEI    = 

.117895E    01 

.1173I4-IE  01 

YT    = 
B   = 

.123598E   01 

.925627E  00 

YE   = 
YEI    = 

.130363E  01 
.129810E  01 

YT    = 
B   = 

.13^85IiE  01 
.113353E  01 

YE   = 
YEI    = 

.1I1.0J133E  01 

.139879E   01 

YT    = 
B   = 

.ll|3773E  01 
.130622E  01 

Table    10    (Concl.) . 


60 


XE  = 
YEE  = 

.9993^-OE  01 
.11^7106E  01 

XT  = 
YDT  = 

.150213E  02 
,14j.1228E-02 

XE  = 
YEE  = 

.225557E  02 
.1526J+9E  01 

XT  = 

YDT  = 

..i^l353E  02 
.LJ.90629E-03 

XE  = 
YEE  = 

.902853E  02 
.155972E  01 

XEN  = 

.902853E  02 

YE  = 
YEI  = 

.ll4.8211|E  01 
.II4.766OE  01 

YT  = 
B  = 

.l50i±32E  01 
.l^oOljE  01 

YE  = 
YEI  = 

.153757E  01 

.153203E  01 

YT  = 
B  = 

,l5^86i|E  01 

.152650E  01 

YE  = 
YEI  = 

.157080E  01 

.156526E  01 

YEN  = 

.157080E  01 
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Table  11.   Asymptotic  10-TRIP. 


EMAX 

= 

,9239l|.lE-02 

XE 

YEE 

= 

.30822liE    00 

.29898J4E  00 

XT 
YDT 

= 

482855E  00 
.810932E  00 

XE 
YEE 

= 

.652152E  00 
.577886E  00 

XT 
YDT 

= 

.826711E  00 

.59ij-Ol8E   00 

XE 
YEE 

= 

.100917E   01 

.78996I4E  00 

XT 
YDT 

= 

.121700E  01 

.I4.O3OI4.7E    00 

XE 
YEE 

= 

.1M|.257E    01 

.96lj.6i4.5E  00 

XT 
YDT 

= 

.171672E  01 

.25331+9E   00 

XE 
YEE 

= 

.202320E  01 

.11117i|E   01 

XT 
YDT 

= 

.2I4.1792E  01 
.11+6063E  00 

XE 
YEE 

= 

.287229E  01 

.123577E    01 

XT 
YDT 

= 

.3i^9535E  01 

.756572E-01 

XE 
YEE 

= 

.i|23720E  01 
.133903E   01 

XT 
YDT 

_ 

.533708E   01 

.339162E-01 

YE   - 
YEI    = 

.30822I|E   00 

.303601JE  00 

YT    = 
B   = 

.I4J4.9838E  00 

.5O2753E-01 

YE   = 
YEI    = 

.587126E  00 
.582506E  00 

YT   = 
B   = 

.690817E  00 

.199736E   00 

YE   = 
YEI    = 

.79920I4E  00 
.79I4.58I4E  00 

YT    = 

B   = 

.882969E  00 
.392I4.59E  00 

YE   = 
YEI    = 

.973885E  00 
.969265E  00 

YT    = 
B   = 

.10li33i4E  01 
.60&J.10E   00 

YE    = 
YEI    = 

.112098E   01 

.111636E  01 

YT   = 
B   = 

.117861|JE   01 

.825470E  00 

YE    = 
YEI   = 

.12l|501E  01 
.12^039E   01 

YT   = 
B   = 

.12921I|E   01 

.102770E  01 

YE   = 

YEI    = 

.13I4-827E  01 

.131^365E   01 

YT   = 
B   = 

.138557E   01 

.120I4.56E  01 

Table    11    (Concl. ) . 
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XE   = 
YEE   = 

.670618E  01 

.1I(.2277E   01 

'    XT    = 
YDT   = 

.897I4-67E  01 

.122632E-01 

XE   = 
YEE   = 

.119931-iE  02 
.HJ.876IE    01 

XT    = 
YDT    = 

.180162E  02 
.307114.1E-02 

XE   = 
YEE   = 

.27010l1E  02 
.l£338l|E   01 

XT   = 
YDT    = 

.514.1082E  02 
.3i4.1I|48E-03 

XE   = 

YEE   = 

.108229E   03 

.156156E  01 

XEN   = 

.108229E   03 

YE   = 

.ll|.3201E   01 

YEI    = 

.11^27392   01 

YT   = 

.1^5983E    01 

B   = 

.13^9 7 7E  01 

YE   = 

.1I4.9685E  01 

YEI    = 

.1[^9223E   01 

YT   = 

.151535E   01 

B   = 

.IJ4.60OIE    01 

YE    = 

.l5i|307E  01 

YEI    = 

.l538i].6E  01 

YT   = 

.155232E  01 

B   = 

.15336I4E  01 

YE    = 

.157080E  01 

YEI    = 

.156618E  01 

YEN   = 

.157080E   01 

Table    12.      Asymptotic   lij.-TRIP. 
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CMAX  = 

.5096i^9E-02 

XE  = 

YEE  = 

.251271E  00 

.2)+6l75E  00 

XT  = 
YDT  = 

.38788L1E  00 

.869222E  00 

XE  = 
YEE  = 

.517036E  00 
.I4.77183E  00 

XT  = 

YDT  = 

.6i|2k99E  00 

.707012E  00 

XE  = 
YEE  = 

.76930!^  00 
.6557J4.2E  00 

XT  = 
YDT  = 

.903179E  00 

.550742E  00 

XE  = 
YEE  = 

.10k25lE  01 
.806210E  00 

XT  = 
YDT  = 

.119669E  01 
.1+11173E  00 

XE  = 
YEE  = 

.136037E  01 

.936902E  00 

XT  = 
YDT  = 

.151|806E  01 
.29l|l+20E  00 

XE  = 

YEE  = 

.175066E  01 
.105181E  01 

XT  = 
YDT  = 

.199056E  01 
.201519E  00 

XE  = 
YEE  = 

.225361E  01 

.115317E  01 

XT  = 
YDT  = 

.257528E  01 
.131026E  00 

XE  = 
YEE  = 

.293389E  01 
.12Lj.230E  01 

XT  = 
YDT  = 

.338813E  01 

.801319E-01 

YE  = 
YEI  = 

.251271E  00 
.2L1.8723E  00 

YT  = 
B  = 

.370019E  00 
.328609E-01 

YE  = 
YEI  = 

.i|82280E  00 
479731E  00 

YT  = 
B  = 

.571083E  00 
.116315E  00 

YE  = 
YEI  = 

.660838E  00 
.658290E  00 

YT  = 

B  = 

.7314.569E  00 

.237150E  00 

YE  = 
YEI  = 

.811306E  00 
.808758E  00 

YT  = 
B  = 

.87^699E  00 

.382653E  00 

YE  = 

YEI  = 

.91+1999E  00 
.9391+51E  00 

YT  = 
B  = 

.997261E  00 

.5I4.1480E  00 

YE  a 
YEI  = 

.105691E  01 

.105^-36E  01 

YT  = 
B  = 

.110525E  01 

.70lj.ll8E  00 

YE  = 
YEI  = 

.115826E  01 

.115571E  01 

YT  = 
B  = 

. 1200klE  01 
.862982E  00 

YE  = 
YEI  = 

.12L^+0E  01 
.12iii|85E  01 

YT  = 

.128380E  01 

B  = 

.101230E  01 

Table    12    (Concl.) . 
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XE   = 
YEE   = 

.390l;12E  01 
.132005E    01 

XT   = 
YDT    = 

.l4.585i4.iE  01 

.li.5i4.OO9E- 01 

XE   = 
YEE   = 

.537717E  01 
.138692E  01 

XT    = 
YDT   = 

.6I4.7893E  01 
.232685E-01 

XE   = 
YEE   = 

.779798E  01 
,li|i325E   01 

XT    = 
YDT    = 

.976967E  01 

.IO368I4E-OI 

XE   = 
YEE   = 

.122312E  02 
.II4.8922E   01 

XT   = 
YDT    = 

.  I6326J4E   02 

.373759E-02 

XE   = 
YEE    = 

.217835E  02 
.152J+92E   01 

XT    = 
YDT   = 

.326896E  02 
.93i|-9l8E-03 

XE   = 
YEE   ■ 

.I4.90I4.53E  02 
.lS&OlJlE    01 

XY   = 
YDT   = 

.981012E  02 
.103898E-03 

XE   = 
YEE   = 

.196208E  03 
.156570E    01 

XEN   = 

.196208E   03 

YE   = 
YEI    = 

.13251I1E  01 
.132259E  01 

YT    = 

B   = 

.135607E  01 

.ll[i789E   01 

YE   = 

YEI    = 

.139202E  01 
.13891j.7E   01 

YT   = 
B   = 

.II4.1766E  01 
.126690E    01 

YE    = 
YEI    = 

.II4I4.835E  01 
.II4I4.58OE   01 

YT   = 

B   = 

.1^6879E   01 

.136750E  01 

YE    = 
YEI    = 

.1i4.9i4.3iE  01 

.1I4.9177E  01 

YT    - 

B   = 

.150962E  01 

.II4I4.86OE   01 

YE   = 

YEI    = 

.153002E  01 
.l527i+7E    01 

YT    = 

B   = 

.1514021E   01 

.150965E  01 

YE   = 
YEI   = 

.155551E   01 

•155296E  01 

YT    = 
B   = 

.156060E  01 
.1550I4.1E  01 

YE   = 
YEI    = 

.157080E  01 
.156825E  01 

YEN   = 

.157080E  01 
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Table  13.   Maximal  errors  for  polygonal  approximations 
with  Chebyshev  norm  and  zero  boundary  error. 

M     :     N     :  EM AX 


1 

1.570796 

2 

0 

0.566912 

3 

1 

0.195652 

k 

2 

0.0991702 

5 

3 

0.0599589 

6 

k 

O.Ol|01706 

7 

5 

0.0287930 

8 

6 

0.0216503 

9 

7 

0.01687258 

10 

8 

0.01351931+ 

11 

9 

0.011075^8 

12 

10 

0.009239^8 

13 

11 

III- 

12 

15 

13 

16 

ik 

O.OO5096I4.9 
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REALIZATION  OF  THE  ARCTANGENT  SPC 


Fundamental  Idea 


This  chapter  conveys  the  basic  idea  for  constructing  an 
SPC,  but  does  not  give  a  hardware  design  of  an  Arctangent  SPC. 
It  will  lead  to  a  better  understanding  of  the  simulation  of  an 
Arctangent  SPC  on  a  G-PC . 

Figure  30  shows  a  schematic  diagram  of  an  SPC  with  upper 
polygonal  approximation.   Next  to  the  input,  which  is  assumed 
to  have  pulse  form,  is  an  encoder  (enc.)  unit  translating  the 
input  signal  x  into  machine  language.   There  the  signal  path 
transmitting  x  splits  up  in  two  paths. 

The  first  path  leads  to  the  decision  network  (dec.  netw.). 
There  the  signal  x  is  compared  at  different  decision  stations 
with  precomputed  signal  levels.   The  number  of  outputs  equals 
the  number  of  pairs  (ni  >b  )  that  are  stored  in  the  nonerasable 
memory  (non-er .mem. ) .   After  having  passed  several  decision 
stages,  the  input  signal  reaches  one  of  the  outputs  of  the  de- 
cision network.   Here  a  driving  pulse  is  generated  and  reads  out 
the  information,  which  thereupon  is  transferred  to  the  arith- 
metic control  (arith.  con.)  unit. 

The  second  path  transmits  x  through  a  delay  (del.)  unit  to 
the  accumulator  (ace),  where  the  arithmetic  operations  speci- 
fied by  the  expression 

y  =  mn  •  x  +  bn 
are  executed.   The  final  result  y  passes  through  a  decoding  (dec. 
unit,  where  y  becomes  the  desired  output  form. 
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non-er. memory 


input 


output 


Pig.  30.   Schematic  diagram  of  an  SPC,  using 
polygonal  approximation. 
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Thus  the  main  parts  of  this  computer  are  the  accumulator, 
the  nonerasable  memory,  and  the  decision  network.   The  preceding 
chapter  included  all  the  numerical  work  used  for  the  realization 
(hardware  design)  or  the  simulation  (software  programming).   The 
following  chapter  takes  up  some  details  about  the  decision  net- 
work. 

The  Decision  Network 

As  argument  x  enters  the  decision  network,  it  is  compared 
to  the  different  precomputed  levels  x  ,,  x  2^  •  •  • >  xeM>  xeN+l* 
Thus  the  appropriate  parameter  pair  (m  ,  b  )  is  read  out  from 
the  immediately  following  nonerasable  memory.   The  comparisons 
can  be  executed  in  many  ways.   Here,  two  extreme  cases,  case 
(a)  and  case  (b),  are  considered. 

Before,  it  should  be  recalled  that  the  maximal  error  argu- 
ments xQ-|_,  xe2,  .  .  .,  xQ^,  xe^+2_  are  the  boundaries  of  the  in- 
tervals I0  b  (0,  xQl),  I-l  =  (x0l,  xe2),  •  •  •  ,  In  =  xeN>  xeN+l> 
X^=    (xeN+l>  xecJ'  ™here  I=I0uIlUl  "  '  w  %  ^  J  oo  =    (°>©°)- 

Case  ( a) --The  Search  Method.   This  method  takes  the  abso- 
lute value  of  the  argument  x  and  saves  its  sign.   Then  it  con- 
tinues with  the  left  most  interval  IQ  =  (0,  xq1)  of  I  =  (0,  OO  ). 
Here  at  the  first  decision  station,  it  is  decided  whether 
x  €  Iq  or  x  6  I-j^  w  I2  ^  •  •  •  ^  Ijj  ^  I  oo  •   In  case  x  £  Iq, 
argument  x  generates  a  pulse  that  is  transmitted  to  the  corre- 
sponding place  in  the  nonerasable  memory  for  read-out  purpose. 
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In  case  x  6l]_  u  I2  u  .  .  .  w  1^  w  1 00  ,    x  is  compared  at  the 
next  decision  station  with  x  2>  where  either  x  e  1^  or 
x  e  Ijv/  .  .  .  u  Ij|  w  I  oo  •   From  here  the  procedure,  as  de- 
scribed for  the  first  decision  station,  starts  afresh. 

If  x  lies  in  the  right  most  interval  1^  ,  x  passes  through 
N+l  decision  levels.   As  it  takes  time  to  execute  these  deci- 
sions, the  processing  time  for  x  £•  1-^  is  much  longer  than  for 
x  £  Iq,  where  only  one  decision  has  to  be  made.   Thus  the  deci- 
sion time  depends  largely  on  the  magnitude  of  the  argument  x. 
This  is  not  an  expedient  feature  if  equal  overall  computation 
time  is  required  regardless  of  the  size  of  the  argument  x.   The 
advantage  of  this  method  is  that  its  decision  process  can  be 
programmed  as  a  loop,  called  "Search  Loop".   Figure  31  shows  a 
flow  chart  of  the  search  method  and  Fig.  32  its  representation 
as  search  loop. 

Case  (b) --The  Tree  Method.   Real-time  control  systems  often 
demand  equal  overall  computation  time  for  their  SPC's,  regard- 
leas  of  the  magnitude  of  the  input  argument.   Th§  computation 
time  (CT)  of  the  ARCTANGENT  SPC  consists  mainly  of  decision  CT 
and  arithmetic  CT.   A  workable  way  to  satisfy  the  demand  of 
equal  overall  CT  is  to  make  decision  CT  and  arithmetic  CT  con- 
stant for  all  possible  arguments  x. 

For  the  arithmetic  unit,  this  is  done  by  not  taking  advan- 
tage of  special  situations,  as  for  x  =  0  or  x  =  00  ,  where  mul- 
tiplication by  1  and  addition  of  0  become  trivial.   These 
special  cases 
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y  -  m0x+b0 


^jy  =  m2x+b2 


y  =  m^x+bo 


Fig.  31.   Flow  chart  of  the  search  method  with 
tree  decision  levels. 


J    1 


J  =  3+1 


Pig.  32.   Flow  chart  for  the  search  loop  method 
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x  =   o   :   y  =  1.00000  •  x  +  0.00000 

x  =  oo  :   y  =  0.00000  •  x  +  1.57080 
are  treated  as  the  others  with  no  such  trivial  values. 

The  decision  CT  can  be  equalized  by  using  the  tree  method. 
Argument  x  passes  through  the  same  number  of  decision  levels, 
regardless  of  the  magnitude  of  x. 

Figure  33    shows  a  decision  tree  with  three  decision  levels. 
This  is  appropriate  if  M  =  2-^  =  8  outputs  are  required,  or  in 
general  M  =  2k  (k  =  1,  2,  3,     ...).   In  case  M  /  2k  (k  =  1,  2, 
3,    ...),  degenerated  decisions  can  be  introduced.   A  degenerated 
decision  element  has  outputs  which  lead  to  the  same  destination. 
(Pig.  3k-) 
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Pig.    3k--      Degenerate    decision. 
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APPENDIX  A 


Calculation  of  the  Argument  with 
Minimum  Curvature 


The  derivatives  of  f(x)  =  arctan(x)  are: 

f'(x)  =  1/(1  +  x2)  (1) 

f"(x)  =  (-2x)/(l  +  x2)2  (2) 

f  »'(x)  =  (6x2  -  2)/(l  +  x2)3  (3) 

The  curvature  of  an  analytic  function  f(x)  is  given  by  the 

equation: 

f"(x) 


curv(x) 


(1  +  (f'(x))2)3/2 
This  equation  derived  with  respect  to  x  yields: 

f'"(x)  •  (1  +  (f'(x))2)3/2  -  3f(x)(f"(x))2 
d(curv(x))  (1  +  {fix))2)1/2 


dx  (1  +  (f'(x))2)3 

The  right  member  of  this  equation  vanishes  at  the  minimum  of 
curvature.   The  inequality 

1  +  (f'(x))2  /  0 
should  be  noticed.   By  taking  into  account  both  remarks,  one 
obtains : 

f"'(x)  •  (1  +  (f'(x))2)  -  3f(x)(f"(x))2  -  0      (ij.) 
Substitution  of  f'(x),  f"(x),  f'''(x)  in  equation  (i|)  by  equa- 
tions (1),  (2),  and  (3)  yields: 
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6x2  -  2          1  1        lpcd 

(1  +  5-0)  -  3 


(1  +  x2)3  *'    (1  +  x2)2      1  +  x2    (1  +  x2)^ 
or 

(6x2  -  2) ( (1  +  x2)2  +  1)  -  12x2  =  0  . 
Further  reduction  yields  the  equation 

x2(6x2  -  2)  (2  +  x2)  -  If.  =  0  . 
The  substitution  x2  =  t  yields 

t(6t  -  2)  (2  +  t)  -    L±  =   0 
or 

6t3  +  lot2  -  l±t   -   l\.   =  0  . 
The  next  chore  is  to  find  the  zero  crossing  of 

h(t)  =  6t3  +  10t2  -  k.t   -  I4.  .  (5)  • 

As  h(0)  =  -I4.  and  h(l)  =  8,  the  zero  crossing  of  (5)  lies  in  the 
interval  (0,  1) .   Further  calculation  yields  an  approximate 
value  for  the  zero  crossing: 

t  =  0.691 
and 

x  =  \l0.691   =  0.831 
is  the  abscissa  value,  where  arctan(x)  has  minimum  curvature. 
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APPENDIX  B 


Sufficient  Conditions  for  Convergence  of 
the  Newton  Process 


Let  "  D "  be  a  general  symbol  for  either  one  of  the  linear 
order  relations  "<  "  or  ".>",  and  let  "□  "  be  a  general  symbol 
for  either  one  of  the  partial  order  relations  "  ^  "  or  "  ^". 
If  "□  "  is  identified  with  "<  ",  then  "□"  is  to  be  identified 
with  "-£  ";  and  if  "O"  is  identified  with  "  ^>  ",  then  "□  " 
is  to  be  identified  with  "^  ". 

Definition  1.   A  function  h(x)  defined  in  an  interval 
( x  ,  x-^)  is  monotonic,  if  for  all  X]_,  X2  £  (xa>  xb)  "k^e  relation 
x-,  □  Xp  implies  h(  x-,  )  □  h( Xp) . 

Definition  2.      A  function  h(x)  defined  in  an  interval 
( xa ,  x-^)  is  strictly  monotonic,  if  for  all  x-|_,  Xp  €  (xQ,  x-^)  the 
relation  x-j_Q  Xp  implies  h(x-]_)  O  h(xp)  .■ 

For  the  next  two  definitions,  the  general  linear  order  re- 
lation "□  "  is  identified  with  the  specific  linear  order  rela- 
tion "  <C  ",  and  the  general  partial  order  relation  "□"  is 
identified  with  the  specific  partial  order  relation  "^  ". 

Definition  3.   A  function  h(x)  defined  in  an  interval 
(x  ,  x,  )  is  increasing  if  for  all  x-,  ,    Xp  £  (x  ,  x-,  )  the  relation 
x-]_  <L  Xp  implies  h(x-j_)  -^h(xp). 

Definition  [j..   A  function  h(x)  defined  in  an  interval 
(xQ,  x-^)  is  strictly  increasing,  if  for  all  x-j_,  Xp  G  (xa>  xt>) 
the  relation  x-j_  <  xp  implies  hCx^)^:  h(xp).   Similarly,  the  dual 
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notions  decreasing  function  and  strictly  decreasing;  function 

are  defined. 

Definition  $.      Let  h(x)  be  a  continuous  function  defined 

in  the  interval  I  =  (xa,  x^)  ,  and  let  Ij  =  (x-]_,  x2)  be  an  open 

interval  such  that  1-^  'SI.   The  points  [x^_,    h(x-^)}  and  (x2,  h(x2)} 

determine  a  secant  line  sv  ^  (x).   The  function  h(x)  is  convex 
x1x2  

in  I  =  (xQ,  x-j^)  ,  if  for  all  intervals  IjCl  and  for  all  x  €  I]_ 
the  inequality  sx  x  (x)  -  h(x)  ^  0  holds. 

Definition  6.  Let  h(x)  be  a  continuous  function  defined  in 
the  closed  interval  I  =  (x  ,  x^)  ,  and  let  I-j_  =  (x-j_,  x2)  be  an 
open  interval  such  that  L  C  I.  The  points  (x-^,  h(x-L)}  and 
\*-2'  ^(x2^)  determine  a  secant  line  sx  x  (x).  The  function 
h(x)  is  strictly  convex  in  I  =  (xa,  xb) ,  if  for  all  intervals 
I  IZ>  I-,  and  for  all  x  6  I-,  the  inequality  sx  x  (x)  -  h(x).>  0 
holds. 

Strict  monotonicity  excludes  horizontal  line  pieces  in 
h(x);  however,  it  does  not  exclude  horizontal  tangent  lines. 
Strict  convexity  excludes  the  equation  of  a  straight  line  as  a 
convex  function.   However,  strict  convexity  is  not  required  for 
the  Newton  process,  as  will  be  seen  later.   All  dual  definitions 
of  terms  like  monotonically  decreasing  function,  concave  func- 
tion, etc.,  as  well  as  the  corresponding  dual  theorems  are 
omitted,  as  both  kinds  are  easily  obtained  by  applying  duality 
rules.   For  the' same  reason,  the  convergence  of  the  Newton  pro- 
cess is  only  proved  for  one  out  of  four  cases  where  h(x)  is 
convex  and  increasing.   In  the  subsequent  theorem,  the  notion 
"continuous  differentiability"  will  appear.   This  is  not  a 
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pleonasm.   The  attribute  "continuous"  excludes  infinite  deriva- 
tives, i.e.,  tangent  lines  that  are  perpendicular  to  the  x-axis. 

Theorem  1.   Let  h(x)  be  a  convex,  continuous,  continuously 
dif f erentiable, and  strictly  monotonic  function  in  the  closed 
interval  (xa,  xb) .   If  h(xa)  <  0  and  h(xb)  >  0,  there  is  a  unique 
real  number  Xq  in  the  interval  (xa,  x^)  for  which  h(xQ)  holds. 
If  the  Newton  process  is  started  at  x  =  x^,  then  it  converges 

tO  Xq  . 

Proof.   Part  A    (Existence  of  Xq) .   Continuity  of  h(x)  and 
the  satisfied  condition  sign  (xa)  ^  sign  (x^)  imply,  by  means 
of  a  theorem  of  Bolzano  that  can  be  found  in  reference  (6,  Bd.  1, 
page  538)  ,'  the  existence  of  an  argument  Xq  in  an  interval 
(xa,  xb)  where  h(xQ)  =  0.   Part  A   Q.E.D. 

Proof.   Part  B  (Uniqueness  of  Xq) .   If  it  is  assumed  that 
h(x)  has  two  zero  crossings,  Xq-^  an^  Xq2,  then  by  Rolle '  s  theorem 
there  exists  an  argument  x^   such  that  Xq->  <.  x^  <  Xq2  and 
h'(x.  )  =  0.   h(*m)  is  the  maximum  of  h(x)  in  interval  (xqt,  xQ2^  ■ 
By  definition  of  such  a  maximum  h(xm),  the  function  h(x)  de- 
creases for  some  x  >  x  .   This  contradicts  that  h(x)  is  a 
strictly  increasing  function.   Hence  h(x)  has  exactly  one  zero. 
Part  B  Q.E.D. 

Lemma  1.   Under  the  assumptions  of  Theorem  1,  h(x) >    0,  if 
x  lies  in  the  interval  (xq,  xb) ,  and  h(x)  >  0  if  x  lies  in  the 
interval  (xa,  xq) . 

Proof.   As  f(xQ)  =  0,  strict  monotonicity  does  not  tolerate 
any  sign  change  inside  the  intervals  (xa,  xq)    and  (xq,  x-5)  . 
Q.E.D. 
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Lemma  2.      Under  the  assumptions  of  Theorem  1,  the  slope 
m(x-,,  Xp)  of  s     ( x)  is  a  monotonically  increasing  function,  if 
s     (x)  is  a  secant  line  through  the  points  {x-,,  h(x-j)  \     and 
{xp,  ^(x2)}  >    where  xi  <   x2*   -^  does  no'fc  matter  whether  x-j_  in- 
creases and  Xp  is  constant  or  Xp  increases  and  x-,  is  constant. 

Proof.   Let  us  consider  the  case  where  Xp  is  constant  and 
x-,  increases  in  the  interval  (xa,  Xp)  .   Lemma  2  is  proved  by 
showing  that  for  two  positions  x-,-,  and  x-,p  of  argument  x-, 
where  x-,-,  <  x-,p,  the  slope  m(x-,-,,  Xp)  is  less  than  the  slope 
m(x-,p,  Xp),  where  x-q  <  x-^p.   The  definition  of  convexity, 
symmetry,  and  monotonicity  is  applied  to  the  secant  line 
sx  x  (x).   Accordingly,  at  argument  x12  holds 

d   =    Slx2Ul2)     "   hUl2)    >  °    '  U) 

The  slope  of  secant  line  sY   v  ( x)  is: 
x12x2 

h(xp)  -  h(x12) 

m(x12,  Xp)  =  .  (2) 

x2  -  x12 

The  slope  of  secant  line  sv   v  (x)  is: 
xllx2 

(h(x2)  -  h(x12))  -  d 

m(xi:L,  x2)  =  .  (3) 

Xp  -  x12 

Since  d  >  0,  and  the  other  differences  are  by  definition  posi- 
tive, it  follows  that 

(k) 

Therefore  in  case  Xp  is  constant,  and  x->    increases.   Lemma  3  is 
valid.   In  the  other  case  where  x-,  is  constant  and  x0    increases, 
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similar  reasoning  leads  to  the  inequality 

m(x1,    x21)  <  mtx-^  x22)  ,  (5) 

where  x-,  -<  x2^  <■  x22  ^x-^.   Q.E.D. 

Lemma  3.   Under  the  assumptions  of  Theorem  1,  the  only 
points  of  intersection  in  the  interval  (xQ,  x-^)  between  the 
function  h(x)  and  the  secant  line  sx  x  (x)  are  the  points 
{x1;,  h(x1)j;   and   {x2,  h(x2)} 

Proof.   It  is  supposed  that  there  is  a  third  intersection 

point  with  an  argument  xin  that  lies  in  the  interval  (xa,  x-j_) 

such  that  xQ  <  x^n  <x-j_.   An  argument  x-^  "begins  to  increase  at 

x-|-j  =  x-  ,  where  the  slope  of  the  secant  line  m(x-,-,,  x-,)  = 

m(xin,  x-j^)  =  m(x1,  x2)  .   As 

ng  to  Lemma 

secant  line  s     (x)  becomes  tangent  line  o"tx^(x)  with  the 

slope  0m(x1).   By  Lemma  2  holds: 

Qm(x1)>  m(x]_,  x2)  .  (1) 

The  existence  of  this  tangent  line  is  guaranteed  by  the  assumed 

continuous  differentiability  of  h(x).   The  proof  goes  on  with  an 

argument  x-,  2  beginning  to  decrease  at  x-^2  =  x2,  where  the  slope 

of  the  secant  line  m(x.  ,  xi?)  =  m(x-,,  x2)  .   x-,2  decreases  until 

it  reaches  x-,  ,  where  the  secant  line  sv   -y  (x)  becomes  tangent 
1'  -x-i  2  2 

line  1tx  (x)  with  slope  ^(x-^).   Prom  Lemma  2  can  be  inferred 
that 

-^(x^)  ■<  m(x-j_,  x2)  .  (2) 

Inequalities  (1)  and  (2)  yield  inequality 

0m(x1)  /  1m(x1)  . 
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This  contradicts  the  previously  assumed  continuous  differentia- 
bility of  h(x) .   Q.E.D. 

From  Lemma  3  follows  immediately  an  interesting  property  of 
continuously  dif f erentiable  functions  h(x).   By  decreasing  x2 
toward  x-,  ,  such  that  after  having  taken  the  limit  x2  — >  x-±, 
secant  line  s„  „  ( x)  becomes  tangent  line  tY  (x).   This  leads 

.A.-]  A.p  -"-~| 

to  Lemma  1+.,    a  modified  form  of  Lemma  3. 

Lemma  l±.      Under  the  assumptions  of  Theorem  1  and  if  tx  ( x) 
is  a  tangent  line  of  h(x)  with  point  of  tangency  at  {xi>    h(x-,  )\  , 
the  only  common  point  of  h(x)  and  tx  (x)  is  the  point  of  tan- 
gency at{x^,  h(x-j_)jin  the  interval  (xQ,  x-^)  such  that  for  all 
x  £  (xQ,  xb)  -   (x-J   holds: 

h(x)  -  t_  (x)  >  0  . 
xl 

Lemma  j?.   Under  the  assumptions  of  Theorem  1,  the  first 
derivative  h'(x)  is  a  monotonically  increasing  function. 

Proof.   Let  us  consider  an  interval  (x-,,  Xp)  that  is  con- 
tained in  interval  (x  ,  x^)  .   The  application  of  Lemma  2  to 
tangent  line  tx  (x)  and  secant  line  sx  x  (x),  and  to  secant  line 
sx  Xo(x)  and  tangent  line  tXp(x)  yields: 


This  implies 


m(tXl)  <  m(sXlX2)  <m(tX2)  . 


m(tXi)<m(tX2) 


As  m(tx  )  =  h'(xo)  and  m( tx  )  =  h'(x2),  the  above  inequality 
yields: 

h'(x1)  <h'(x2)  . 
In  other  words,  x^_  <  x2  implies  h' (  X]_)  -<  h' (  x2)  .   Q.E.D. 
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Lemma  6.   Under  the  assumptions  of  Theorem  1,  zero  xQ  of 
the  function  h(x)  is  a  zero  crossing. 

Proof.   By  Lemma  $,    x,  the  derivative  of  h(x)  at  xq,  is 
different  from  zero.   Q.E.D. 

Lemma  "]_.      Under  the  assumptions  of  Theorem  1,  the  only  zero 

derivative  h'(x)  =  0  occurs  at  argument  x  ,  and  the  inequality 

h'(x)  >  0  holds  in  the  left  half  interval  (xa,  xb) . 

Proof.   Let  us  consider  the  secant  line  s„  v    ( x) ,  where 
xaxl 

x  <C  x-i  *<Xi^.   The  slope  of  this  secant  line  is: 

h(xa)  -  h(xx) 

m(xa,  x1)  =  

xa  -  Xl 

The  strictly  increase  of  h(x)  implies  that  the  above  quotient 

is  positive  for  all  x^'s  in  left  open  interval  (xQ,  xfe)  .   As  x^_ 

tends  to  xa,  one  obtains 

h(xa)  -  h(x-j_) 

lim    =  h'(xa)  —  ° 

xl^  xa      xa  "  xl 

The  limit  has  the  possibility  to  take  to  the  value  zero.   This 
proves  the  first  assertion.   As  h'(x)  is  strictly  increasing 
by  Lemma  5  in  closed  interval  (x  ,  xO  and  as  h'(x  )  =t  0,  h(x) 
is  positive  throughout  interval  (x  ,  xfe) .   Q.E.D. 

Lemma  8.   Under  the  assumptions'  of  Theorem  1,  h'(x)  is 
never  infinite  throughout  (xa,  x-^)  . 

Proof.   This  follows  immediately  from  continuous  differ- 
entiability, which  excludes  h'(x)  =  c*o  . 

Lemma  9.   If  argument  x-^  lies  in  the  left  open  interval 
(xq,  x^) ,  where  xQ  is  the  zero  crossing  of  h(x),  then  the 
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tangent  line  tx  (x)  intersects  the  x-axis  at  X2  such  that 
Xq  <  x2  <  xx. 

Proof.   By  Lemma  Ij.  holds  h(x)  -  tx  ( x)  >  0.   Therefore 
h(x2)  -  tx  (x2)  >  0.   Here  Part  A  and  Part  B  of  Theorem  1  are 
applied  to  interval  (x  ,  Xp)  .   According  to  these  proved  asser- 
tions, there  is  a  unique  zero  crossing  Xq  in  the  interval 
(x  ,  x-u)  .   Continuous  differentiability  of  h(x)  prohibits 
xq  =  x2  by  Lemma  8.   As  h(xp)  •>  0  and  h'(x)  >  0,  the  relation 
Xq  <  xp  holds.   The  same  reasoning  applied  to  argument  x-j_  leads 
to  the  inequality  Xp  <:  x^_.   Q.E.D. 

Proof.   Part  C  (Convergence  of  Newton  Process) .   Let  us 

start  the  procedure  known  as  Newton  process  at  x  =  x-j-,  =  X]_. 

The  intersection  point  Xp  of  the  tangent  line  tx  (x)  and  the 

x-axis  are  given  by 

h(xx) 

x?  =  Xl . 

h'(xx) 

Xp  lies  in  the  open  interval  xq  <  x2  <  x^  by  Lemma  9.   Inequality 
Xq  <  Xp  and  strict  increase  of  h(x)  imply  h(xp)  >  0.   Therefore 
this  procedure  can  be  applied  again  to  argument  xp .   This  itera- 
tive procedure,  called  a  Newton  Process,  yields  the  recurrence 
relation 

h(Xl) 

xn+l  =  *n  "  — "      (n  =  1,  2,  3,  .  .  .  ) 

h' (xx) 

The  repetitive  execution  of  this  relation  yields  the  chain 
inequality 

Xq  <  .  .  .  x^  <  xp  <  x1 
where  x-j_  is  the  starting  point  of  the  Newton  process.   By 
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choosing  a  sufficiently  high  value  for  n,  the  above  chain  in- 
equality enables  to  satisfy  the  relation 

*n  "  x0  <  £• 
whichever  small  positive  value  e  has.   In  other  words, 

llXD   ^  =  XQ  . 
n  — ►  <*> 

This  completes  the  proof  of  Part  C  of  Theorem  1. 

Definition  7.   A  Newton  process  for  a  function  f(x)  de- 
fined in  the  interval  (xg,  x-^)  and  with  a  zero  crossing  at  xq 
such  that  xQ<r(x  ,  x,  )  ,  is  strictly  convergent,  if  it  converges 
to  Xq  for  any  starting  point  x,  €  (eo>  xv^  * 

Assumptions  of  Theorem  1  are  sufficient  for  strict  con- 
vergence of  the  Newton  process  in  interval  (x  ,  x.) .   It  should 
be  observed  that  the  above  described  Newton  process  does  not 
make  use  of  any  assumption  of  Theorem  1  for  the  interval 
(x  ,  Xq) .   Therefore  Theorem  1  is  restated  with  slightly  weaker 
assumptions. 

Theorem  2.      If  a  function  h(x)  has  a  zero  and  a  positive 
derivative  at  Xq  and  if  h(x)  is  continuous,  convex,  strictly  in- 
creasing, and  continuously  diff erentiable  in  the  interval 
(xq,  x^) ,  then  the  Newton  process  is  strictly  convergent  for 
interval  (x  ,  xO . 

It  should  be  noted  that  in  the  numerical  computation  of  a 
Newton  process,  assumptions  stated  above  might  be  satisfied  and 
yet  the  calculation  might  become  unstable.   The  reason  for  this 
instability  is  that  the  above  theory  does  not  take  account  of 
the  finiteness  of  registers  and  computation  time. 


APPENDIX  C 

FORTRAN  Program  of  Asymptotic  N-Trip  with  Newton 

Process  in  E-Start,  T-,  and  E-Algorithm  and 

with  Unilateral  Dit-by-dit  Method  in 

Pi  Half  Algorithm 


MCNSS      JOB   ASYMPTOTIC   N  TRIP  ( NP )  UMILAT  APPR 

MCNSS      CCMT  10  MINUTES*  10  PAGES 

MCNSS      ASGN  MJB*12 

MCNSS      ASGN  MGC16 

MCNSS      MODE  GO* TEST 

MCNSS      EXEO  FORTRAN, ,, 10,03, , *ASYMPNTRIP 

800  FORMAT ( 13) 

801  FORMAT ( 1H1»30H  ASYMPTOTIC   N  TRIP(NP)   N  =   ,13) 

803  FORMAT* 1HL*10H    EMAX  =  iE20.14>10H    DELT  =  >E20.14) 
805  FCRMAT( 1HT,10H   DELTA  =  .E20.14) 

802  FCRMATdOX  ,E20.14,10X  ,E20.14) 

804  FORMAT ( 10X  *E20.14) 

820  FORMAT ( 1HK,10H 

821  FORMAT ( 1HT,10H 

830  FORMAT ( 1H  ,10H 

831  FORMAT( 1HS.10H 
855  FORMAT! 1HT*10H 

860  FORMAT! 1HL,10H 

861  FORMAT! 1HL»10H 
READ  (1,800) 
WRITE  (3,801) 
READ  (1,802) 
WRITE  (3*803) 
READ  (1,804) 
WRITE  (3,805) 

A  =  1.0 
1  WRITE  (3*860)  E 
CALLFLTPT(EMAX) 
XT  =  0.0 


XE  = 

»E2G< 

,14,10H 

YE  = 

♦E20.14) 

YFE  = 

»E20, 

,14»10H 

YEI  = 

♦E20.14) 

1      XT  = 

»E20, 

.  14,10H 

YT  = 

♦E20.14) 

YDT  = 

»E20, 

,14»10H 

D  = 

•E20.14) 

XEN  = 

>E20« 

,14»10H 

YEN  = 

.E20.14) 

l    EMAX  = 

,E20, 

.14) 

i  EMAXEN  = 

N 

»E20. 

,14) 

N 

EMAX,  DELT 

EMAX*  DELT 

DELTA 

DELTA 

YT  = 

0.0 

YDT  = 

1.0 

B  = 

0.0 

XE  = 

A 

YEE  = 

ATAN!XE) 

, 

EPSXE  = 

(XE-YEE) 

-  EMAX 

IF(EPSXE) 

XE  = 

XE  +  A 

GO  TC    2 

YDE  = 

1 .  /  (  1 .  + 

XE**2. ) 

IF(ABS(EPSXE)  -  DELTA) 

XE  = 

XE  -  EPSXE/(1.0-YDE) 

YEE  = 

ATAN(XE) 

EPSXE  = 

(XE-YEE) 

-  EMAX 

GO  TO    4 

YE  = 

XE 

YEI  = 

•5*<YE+YEE ) 

WRITE(3»820)    XE* 

YE 

WRITE(3*821)   YEE* 

YEI 

12, 
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DC   91  J=1.N 

XT  =  XE  +  A 

22  YDT  -  l./( 1  .  +  XT**2. ) 

YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT*(XT-XE) 
IF(EPSXT)  23*   25,   25 

23  XT  =  XT  +  A 
GO  TO  22 

25  IF(ABSCEPSXT)  -  DELTA)  32,   32,'   26 

26  XT  =  XT  -  EPSXT/(2.*XT*(XT-XE )*( YDT**2. ) ) 
YDT  =  l./( 1.  +  XT**2. ) 

YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT*(XT-XE) 
GC  TO   25 
32        B  =  YT  -  YDT*XT 
WRITE(3,830)    XT,    YT 
WRITE(3,831)   YDT,     B 
XE  =  XT  +  A 

42  YEE  =  ATAN(XE) 

EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
IF(EPS) I)  43,   44,   44 

43  XE  =  XE  +  A 
GC  TC   42 

44  YDE  =  l./(  1.  +  XE**2. ) 
IF(ABS(EPSXE)  -  DELTA)  52,   52,   46 

46       XE  =  XE  -  EPSXE/(YDT-YDE  ) 
YEE  =  ATAN(XE) 
EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
GC  TC   44 
52       YE  =  YDT*XE  +  B 
YEI  =  .5*(YE+YEE  ) 
WRITE(3,820)    XE,    YE 
WRITE(3,821)   YEE,   YEI 

DEL  =  YE  -  1.570  796  326  794  9 
IF(DEL)  91,   91,   63 

91  CCNTINUE 

WRITE  (3,855)  XE,  YE 

IF(ABS(DEL)  -  DELTA)  192,  192,   61 

61  IF(DEL)  62,   63,   63 

62  EMAX  =  EMAX  +  DELT 
GC  TC    1 

63  EMAX  =  EMAX  -  DELT 
DELT  =  .1*DELT 
EMAX  =  EMAX  +  DELT 

GC  TC    1 
192  WRITE  (3,861)  EMAX 
END 
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1 

MON$$ 

EXEQ 

AUTOCODER,, 

»NC 

^AC 

2 

HEADR***FLCATING 

PC 

INT 

TYPE 

SUBROUTINE*** 

3 

TITLEFLTPT 

4 

SBR 

X12 

5 

BCE 

END, /MCI/ 

6 

MLCA 

4+X13»*+6 

7 

MLC 

0,WORK 

7] 

MLC 

8 

BZN 

*-ll,/CTB/ 

9 

BXPA 

/CNC/ 

10 

DCW 

MESSAGE 

12 

END 

B 

5  +  X13 

13 

MESSAGE 

DCW 

=  10 

14 

WORK 

DCW 

=  2 

16 

DCW 

-)  - 

17 

MON$$ 

END 

EXEQ 

CALL 

LINKLOAD 
ASYMPNTRIP 

MONSS 

EXEQ 

ASYMPNTRIP, 

^ji; 

MCN$$ 

JOB 

ACTSS    BEAT  A 

.  GIMMEL 

EE       0 

0154U40408 

The  flow  chart  of  this  program  is  given  in  Pig.  29.   The 
name  of  the  above  AUTOCODER  subprogram  is  PLTPT.   Each  time  the 
statement  FLTPT(EMAX)  is  encountered  in  the  main  program,  the 
actual  value  of  EMAX  appears  on  the  console  typewriter  of  the 
IBM  llj.10  Data  Processing  System. 
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APPENDIX  D 

FORTRAN  Program  of  Asymptotic  N-Trip  with  Newton 

Process  in  E-Start,  T-,  and  E-Algorithm  and 

with  Oscillatory  Dit-by-dit  Method  in 

Pi  Half  Algorithm 
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MCNSS      JOB   ASYMPTOTIC   N  TRIP  (  NP  )  CSCILL  APPR 

MCNSS      CCMT  10  MINUTES*  10  PAGES 

MCNSS      ASGN  MJB.12 

MCNSS      ASGN  MGC16 

MCNSS      MODE  GO.TEST 

MCNSS      EXEQ  FCF TRAN » , , 10 .03 > , » ASYMPNTR I P 

800  FORMAT ( 13 ) 

801  FORMAT ( 1H1.30H  ASYMPTOTIC   N  TRIP(NP)   N  =   ,13) 

803  FORMAT ( 1HL.10H    EMAX  =  .E20.14.10H    DELT  =  ,E20.14) 
805  FORMAT( 1HT.10H   DELTA  =  ,E20.14) 

802  FORMATdOX  .E20.14.10X  .E20.14) 

804  FCRMAT( 10X  ,E20.14) 

820  FORMAT( 1HK.10H      XE  =  .E20.14.10H      YE  =  .E20.14) 

821  FCRMAT( 1HT.10H     YEE  =  .E20.14.10H     YEI  =  .E20.14) 

830  FORMATtlH  ,10H      XT  =  .E2C.14.10H      YT  =  .E20.14) 

831  FORMATS 1HS.10H     YDT  =  .E20.14.10H       B  =  .E20.14) 
855  FORMAT( 1HT.10H     XEN  =  .E20.14.10H     YEN  =  .E20.14) 

860  FORMAT* 1HL.10H    EMAX  =  , £20,14) 

861  FORMAT* 1HL.10H  EMAXEN  =  .E20.14) 
READ   (1.800)  N 

WRITE  (3.801)  N 
READ   (1,802)  EMAX,  DELT 
WRITE  (3,803)  EMAX,  DELT 
READ   (1,804)  DELTA 
WRITE  (3,805)  DELTA 
A  =  1.0 

1  WRITE  (3,860)  EMAX 

XT  =  0.0 

YT  =  0.0 

YDT  =  1.0 

B  =  0.0 

XE  =  A 

2  YEE  =  ATAN(XE) 

EPSXE  =  (XE-YEE)  -  EMAX 
IF(EPSXE)  3,    4,    4 

3  XE  =  XE  +  A 
GO  TO    2 

4  YDF  =  l./C  1«  +  XE#*2.  ) 
IF(ABS(EPSXE)  -  DELTA)  12,   12,    6 

6       XE  =  XE  -  EPSXE/( 1.0-YDE) 
YEE  =  ATAN(XE) 
EPSXE  =  (XE-YEE)  -  EMAX 
GO  TO    4 
12      \:  =  XE 

YEI  =  a5*(YE*fEE) 
WRITE(3.820)    XE,    YE 
WRITE(3,821)   YEE.   YEI 
DO   91  J=1,N 

XT  =  XE  +  A 
22      YDT  =  l./(l.  +  XT**2.) 
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YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT*(XT~XE) 
IF(EPSXT)  23,   25,   25 

23       XT  =  XT  +  A 
GO  TC   22 

25  IF(ABSIEPSXT)  -  DELTA)  32,   32,   26 

26  XT  =  XT  -  EPSXT/(2.*XT*(XT-XE )*(YDT**2.) 
YDT  =  l./< 1.  +  XT**2. ) 

YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT* ( XT-XE ) 
GO  TC   25 
32        B  =  YT  -  YDT*XT 
WRITE(3,830)    XT,    YT 
WRITE(3,831)   YDT,     B 
XE  =  XT  +  A 

42  YEE  =  ATAN(XE) 

EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
IF(EPSXE)  43,   44,   44 

43  XE  =  XE  +  A 
GC  TC   42 

44  YDE  =  l./(  1.  +  XE**2.) 
IFtABS(EPSXE)  -  DELTA)  52,   52,   46 

46       XE  =  XE  -  EPSXE/(YDT-YDE) 
YEE  =  ATAN(XE) 
EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
GC  TC   44 
52       YE  =  YDT*XE  +  B 
YEI  =  .5*<YE+YEE) 
WRITE(3,820)    XE,    YE 
WRITE(3,821)   YEE,   YEI 

DEL  =  YE  -  1.570  796  326  794  9 
IF(DEL)  91,   91,   63 

91  CONTINUE 

WRITE  (3,855)  XE,  YE 

IF(ABS(DEL)  -  DELTA)  192,  192,   61 

61  IF(DEL)  62,   63,   63 

62  EMAX  =  EMAX  +  DELT 
GC  TC    1 

63  DELT  =  0.1  *  DELT 
EMAX  =  EMAX  -  DELT 

101  WRITE  (3,860)  EMAX 
CALLFLTPT(EMAX) 

XT  =  0.0 

YT  =  0.0 

YDT  =  1.0 

B  =  0.0 

XE  =  A 

102  YEE  =  ATAN(XE) 
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EPSXE  =  (XE-YEE)  -  EMAX 
IF(EPSXE)  103,  104,  104 

103  XE  =  XE  +  A 
GO  TO  102 

104  YDE  =  1«/(1*+XE**2«) 

IF(ABS(EPSXE)  -  DELTA)  112,  112,  106 

106       XE  =  XE  -  EPSXE/U.O-YDE) 
YEE  =  ATAN(XE) 
EPSXE  =  (XE-YEE)  -  EMAX 
GO  TO  104 
112       YE  =  XE 

YEI  =  .5*(YE+YEE) 
WRITE(3,820)    XE*    YE 
WRITE(3*821)   YEE,   YEI 
DO  191  J=1,N 

XT  =  XE  +  A 

122  YDT  =  !•/( 1.  +  XT**2. ) 

YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT*(XT-XE) 
IF(EPSXT)  123*  125,  125 

123  XT  =  XT  +  A 
GO  TO  122 

125  IF( ABS(EPSXT)  -  DELTA)  132*  132,  126 

126  XT  =  XT  -  EPSXT/(2.*XT*(XT-XE)*(YDT#*20 > 
YDT  =  l./(l«  +  XT**2. ) 

YT  =  ATAN(XT) 
EPSXT  =  (YT-YE)  -  YDT*(XT~XE) 
GO  TO  125 
132        B  =  YT  -  YDT*XT 

WRITE(3,830)    XT,    YT 
WRITE(3,831)   YDT,     B 
XE  =  XT  +  A 

142  YEE  =  ATAN(XE) 

EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
IF(EPSXE)  143,  143,  144 

143  XE  =  XE  +  A 
GO  TO  142 

144  YDE  =  l./(l.  +  XE**2«) 
IF(ABS(EPSXE)  -  DELTA)  152,  152,'  146 

146       XE  =  XE  -  EPSXE/(YDT-YDE) 
YEE  =  ATAN(XE) 
EPSXE  =  ((YDT*XE  +  B)  -  YEE)  -  EMAX 
GO  TO  144 
152       YE  =  YDT*XE  +  B 
YEI  =  .5*(YE+YEE) 
WRITE(3,820)    XE,    YE   ' 
WRITE(3,821)   YEE,   YEI 

DELN  =  YE  -  1.570  796  326  794  9 
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IF(DELN)  191,     191,     161 

191  CONTINUE 

WRITE' (3,853)  XE,  YE 

161  IF(DELN  *  OtL)  162,  162,  163 

162  DELT  =  .1  *  DELT 

163  DEL  =  DELN 

IF(ABSCDEL)  -  DELTA)  192,  192,  164 

164  IF(DEL)  165,  166,  166 
16  5     EMAX  =  EM AX  +  DELT 

GO  TO  101 
166     EMAX  =  EMAX  -  DELT 
GO  TO  101 

192  WRITE  (3,861)  EMAX 
END 
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APPENDIX  E 
Simulation  of  an  Arctangent  SPC  on  IBM  1620 

The  computation  program  that  is  going  to  be  described  is 
written  in  SPS  -  language  for  electronic  digital  computer  IBM 
1620.   The  SPS  -  processor  program  that  translates  SPS  - 
language  into  the  machine  language  of  IBM  1620  is  an  assembler 
program  such  that  arithmetic,  branch,  and  input-output  SPS 
statements  correspond  in  a  unique  manner  to  IBM  1620 ' s  machine 
instructions.   This  one-to-one  correspondence  and  SPS's  mnemonic 
features  predestinate  this  symbolic  program  language  to  be  used 
for  simulating  an  arctangent  SPC  on  IBM  1620. 

The  subsequent  program  consists  of  two  main  programs  A  and 
B.   Program  A  computes  the  arctangent  (ATAW)  function  by  means 
of  an  upper  polygonal  approximation,  whereas  program  B  computes 
the  ATAN  function  by  using  Chebyshev  power  series.   The  upper 
polygonal  approximation  of  program  A  employs  l6-line  segments 
whose  upper  boundaries  XE,  slopes  YDT,  and  y-intercepts  B  were 
calculated  in  the  Asymptotic  Ii4.-Trip  (cf  Table  12).   In  this 
program  the  triple  (XE,  YDT,  B)  is  denoted  by  (XHIGH,  MSLOPE, 
BINTCP) .   The  Chebyshev  power  series  approximation  of  program  B 
calculates  arctan(x)  for  the  interval  (0,1)  in  a  straight- 
forward manner,  whereas  for  (1,  00  )  arctan(l/x)  is  computed 
first,  then  the  transformation  arctan(x)  =  7t/2  -  arctan( l/x) is 
used.   An  exhaustive  description  of  the  IBM  1620  machine  in- 
structions and  the  SPS  program  language  is  given  by  (8). 
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0 

o2l4u 

DCRG*-3 

02150C 

TFM 

NEFRTN, 00, 10 

02160 

ba 

02170 

DCRG*-9 

02180NEFRTNDS 

10 

02190E14N/ATF 

AEFF  TN, AFIELD 

02200 

BNF 

E,NEFRTN-2, , MANTISSA  SIGN 

O3OIO 

TDM 

AEFRTN-27,2,11 

0?020 

CF 

NEFRTN-2 

O3030E 

CF 

NEFRTN-9 

03040 

TD 

AEFRTN-22, NEFRTN-9 

O3050 

TD 

AEFRTN-20,NEFRTN-8 

O3O6O 

TD 

AEFRTN-18 ,NEFRTN-7 

O3070 

TD 

AEFRTN-16»NEFRTN-6 

O3O8O 

JO 

AEFRTN-14,NLFRTN-5 

03090 

TD 

AEFRTN-12»NEFRTN-4 

03100 

TD 

AEFRTN-10.NEFRTN-3 

0-2  110 

TD   AEFRTN-8,NEFRTN-2 

0  7l?0 

PNF  F, NEFRTN 

0-5130 

TDM  AEFRTN-5,2 

03140 

CF   NEFRTN 

03150F 

CF   NEFRTN-1 

0?160 

TO   AEFRTN-2, NEFRTN-1 

0?170 

TO   AEFRTN»NEFRTN 

0  3  1 8  0 

BB 

0^190 

D0RG*-9 

07200ERRCR  DAC  14, ERROR  IN  DATA- 

0  4010AFIELDDC   28 »000 370 70707070 707 07045 107070 

04020STARTSRCTY 

04030       WATYCMNTA 

04040G      H 

04050       BC1  PROGMB 

O4O6OPROGMARCTY 


04070 

WATYCMNTB 

04080 

TFM 

PRCGMA+6»16»10 

04090 

TFM 

AA+6»XHIGH>7 

0  4  1  0  0 

TFM 

AB+6,MSL0PE»7 

04110 

TFM 

AC+6,3INTCP»7 

O4 120 

TFM 

ADRESA+6»AD,7 

04130AD 

RACDINPUT 

04140 

SF 

INPUT-1 

04150 

SF 

INPUT+27 

0  4  1  6  0 

SF 

INPUT+55 

04170 

BT 

E14A/N, INPUT+26 

O4I8OAA 

TFL 

0, NEFRTN 

04190 

BT 

E14A/N, INPUT+54 

04200AB 

TFL 

0, NEFRTN 

O5OIO 

BT 

El4A/N,INPUT+82 

05O20AC 

TFL 

0, NEFRTN 

0*030 

AM 

AA+6, 10,10 

O5O4O 

AM 

AB+6, 10,10 

0<^050 

AM 

AC+6,10,10 

0*060 

SM 

PRCGMA  +  6,01  ,10 

0^070 

BH 

AD 

0  5  0  8  0 

RCTY 

0*090 

WATYCMNTCA 

05100H 

RCTY 

0  5 1 1 0 

RNTYSTARTS+2» » » TIMES  LC 

OP  NUMBER 

05120 

BC4 

H' 

05130 

SF 

STARTS+2 

0  5 1.40 

RCTY 

OS  150 

WATYCMNTC 

05160 

TFM 

ADRESA+6,SWA,7 

05170J 

RACDTNPUT 

05180 

SF 

INPUT-1 

0*190 

BT 

E14A/N, INPUT+26 

0S200 

TFL 

XNUMBR, NEFRTN 

100 


06010  TOM  H  +  7,0 

0f020  BNF  AF,XNUMBR-2 

06030  SF   H  +  7 

06040  CF   XNUMBR-2 

0A050AF  TF   H+6,STARTS+6 

06051  H 

060601  TFL  G+ll.XNUMBR 

06070  FSUBG+11»XHIGH+150 

06080  BNH  XA 

06090  BTFLE14I^  /A,XNUMBR 

06100  TF   INPUT+26,AEFRTN 

06110  TR   INPUT+27,CMNTD-1 

0*120  WACDINPUT 

06130  RCTY 

06140  WAT Y INPUT 

06150  B    SWA 

06  160  DCRG*-3 

06170XA  TFL  G+11»XNUMBR 

06180  FSUBG+11»XHIGH+   70 

06190  Bri   XI 

06200  TFL  G+11»XNUMBR 

O7OIO  FSUBG+11,XHIGH+   30 

O7O2O  BH   XE 

07030  TFL  G+11»XNUMRR 

07040  FSUBG+11 ♦XHIGH+   10 

07050  BH   XC 

07060  TFL  G+11»XNUMBR 

07070  FSUBG+ll.XHIGH 

O7O8O  BH   XB 

O7O9O  TFM  AH+11,MSLCPE 

071C0  TFM  AI  +  11. BINTCP 

07110  B    AK 

07120  DCRG*-3 

07130XB  TFM  AH+ll.MSLCPE+10 

07140  TFM  AI+11 .BINTCP+lO 

07150  B    AK 

07160  D0RG*-3 

07170XC  TFL  G+ll.XNUMBR 

07180  FSUBG+11. XHIGH+  20 

071<?C  BH   XD 

07200  TFM  AH+11  .MSLCPE+20 

OpOlO  TFM  AI+11, BINTCP+20 

08020  B    AK 

08  03  0  DCRG*-3 

0R040XD  TFM  AH+ 1 1 »MSLCP E+30 

08050  TFM  AI  +  11  .BINTCP  +  30 

0R060  R    AK 

0R070  DCRG*-3 

0P080XE  TFL  G+ll.XNUMRR 

08090  FSUBG+11 .XHIGH+  50 
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0  8  1 0  0 
0  fi  1 1 0 
0  R 1  2  0 

0«130 
08140 

Or  150 

Op  160 

0817  0 

08180XF 

0P19C 

0B200 

0Q010 

0O02CXG 

000  30 

Oo040 

Oq050 

O0O6O 

0o070 

O9O8O 

09090XH 

0Q100 

09110 

0Q120 

0Q130XI 

0ol40 

0o]50 

0Q160 

0Q170 

O0I8C 

0  9  19  0 

09200 

100 10 

10020 

10030 

10040 

10050 

10060XJ 

10070 

10080 

10090 

10100XK 

10110 

10120 

10130 

10  140 

10150  ' 

1016  0 

10170XL 

10180 

10190 


BH   XG 
TFL  G+1  1 
FSUBG+11 
BH   XF 
TFM  AH+1 
TFM  AI  +  1 
B    AK 
DCRG*-3 
TFM  AH+1 
TFM  AI+1 
B    AK 
DCRG*-3 
TFL  G+ll 
FSUBG+11 
BH   XH 
TFM  AH+1 
TFM  AI+1 
B    AK 
DCRG*-3 
TFM  AH+1 
TFM  AI+1 
B    AK 
D0RG*-3 
TFL  G+ll 
FSUBG+11 
BH   XM 
TFL  G+ll 
FSUBG+11 
BH   XK 
TFL  G+ll 
FSUBG+11 
BH   XJ 
TFM  AH+1 
TFM  AI+1 
B    AK 
DCRG*-3 
TFM  AH+1 
TFM  AI+1 
B    AK 
DCRG*-3 
TFL  G+ll 
FSUBG+11 
BH   XL 
TFM  AH+1 
TFM  AI+1 
B    AK 
DCRG*-3 
TFM  AH+1 
TFM  AI+1 
B    AK 


XNUMRR 
XHIGH+  40 


1 ♦MSLCPE+40 
l»BINTCP+40 


l»MSLCPE+50 
1  ♦BINTCP+50 


♦XNUMBR 

»XHIGH+  60 


1 ♦MSLCPE+60 
l»3INTCP+60 


l»MSLCPE+70 
l»BINTCP+70 


♦XNUMBR 
♦XHIGH+  110 


♦XNUMRR 
♦XHIGH+  90 


iXNUMBR 
,XHIGH+  80 


1 ♦MSLCPE+80 
l»SINTCP+80 


l.MSLCPE+90 
l»BINTCP+90 


♦XNUMBR 
♦XHIGH+  100 


1»MSLCPE+100 
l.BINTCP+100 


1  ♦MSLCPE+110 
1  ♦BINTCP  +  IIO 
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ln?00       DCRG*-3 
11010XM     TFL  G+l] ,  XNUMBR 
n02P  FSUBG+1  1  ,XHIGH+  130 

1103C        BH   XC 

nO'iO  TFL    G+ll  ,  XNUMBR 

H050  FSUBG+ll,XHIGh+    120 

11060  BH       XN 

11070       TFM  AH+11 »MSLCPE+120 

11080       TFM  AI  +  H  ,BI.\TCP  +  120 

H090       B    AK 

11 100       DCRG*-3 

llllOXN  TFM    AH+11, MSLOPE+130 

11  120  TFM    AI  +  11,BINTCP+130 

11  130  B         AK 

11140  DCRG*-3 

U  15CXC     TFL  G+l  1  .XNUMRR 

I1  160  FSUBG+11 ,XHIGH+  140 

11170  BH   XP 

11.18C  TFM  AH+11  ,MSLOPE+140 

1U90  TFM  AI  +  H,BINTCP+140 

11200  B    AK 

12010  DCRG*-3 

1?020XP  TFM  AH+11 tMSLCPE+150 

1?03C  TFM  AI+ll.BINTCP+150 

12040  R    AK 

12050  DCRG*-3 

12060AK  TFL  YNUMBR , XNUMBR 

12070AH  FMULYNUMBR»0 

12080AI  FADDYNUMBR,0 

12  090  SM   H+6,01,10 

12  100  BH   I 

12  101  H 

12110  BNF  AL,H+7 

12  120  SF   XNUN 3R-2 

12  130  SF       YNUViflR-2 

12140AL  TR       I\'PUT-1,ZERC-1 

12150  TR   INPUT+79,ZERC 

1?1^C  BTFLF14N/A, XNUMBR 

T170  TF   INPUT  +  26»AEFRTM 

12180  BTFLE14N/A, YNUMBR 

12190  TF   INPUT+56*AEFRTN 

17  200  WACD INPUT 

13010SWA  BNC1J 

13020PRCGMBRCTY 

13030  WATYCMNTCA 

1^040TCERR  RCTY 

13050  RNTYSTARTS+2»», TIMES  LOOP  NUMBER 

1^060  PC4  TCERR 

1^070  SF   STARTS+2 

1 ^080  RCTY 
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nOQG       WATYCMNTC 

13100       TFM  ADRESA+6,SWB,7 

1 -all  OAM     RACDINPUT 

1**12C        SF   INPUT-1 

13130       BT   E14A/N, INPUT+26 

13140       TFL  XNUMBR»NEFRTN 

13150  TF   H+6.STARTS+6 

13151  H 

13160K      FATNYNUMBR»XNUM8R 

13170       SM   H+6»01,10 

1?180       BH   K 

1^181        H 

1^190       TR   INPUT-1,ZERG-1 

1^200       TR   INPUT+79»ZERC 

14010       BTFLE14N/A,XNUM8R 

14020       TF   INPUT+26, AEFRTN 

14030       BTFLE14N/A»YNUMBR 

I4O4O       TF   INPUT+56>AEFRTN 

14050       WACDINPUT 

I4O6OSWB    BC1  AM 

14070       B    PRCGMA 

I4O8O       D0RG*-3 

14090CMNTCADAC  33»ENTER  NUMBER  OF  FUNCTION  PASSES. - 

14100INPUT  DAS  80 

14110       DS   2 

14120CMNTA  DAC  48»SW1  ON  FOR  ARCTAN»CFF  FOR  FUNCTION  APROX I  MAT  I ON- 

14130CMNTB  DAC  33, ENTER  FUNCTION  APROXIMATION  DATA- 

14140CMNTC  DAC  8, ENTER  X- 

14150XHIGH  DSB  10,16 

14160MSL0PEDSB  10,16 

1417CBINTCPDSB  10,16 

14180XNUMBRDS   10 

14  190YNUMBRDS   10 

1420CCMN*D  DAC  34,  =X  IS  NOT  IN  RANGE  OF  FUNCTION.  - 

15010ZERO   DAC  41, 

1*020NCISE  DAC  1,0 

15030       DENDSTARTS 
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PRECCMPUTED 


rUNCTICN   CONSTANTS 


+.2512 
+.5170 
+.7693 
+.1042 
+.1360 
+.17  50 
+.2253 
+.2933 
+.3904 
+.5377 
+.7797 
+.1 223 
+.2178 
+.4904 
+.1962 
+.9999 


7153F+00+, 
3588E+00+, 
0414E+00+, 
51'  0E+01+. 
3676E+G1+, 
6624E+0I+, 
6126E+G1+, 
893lE+ol+, 
12] 3E+01+. 
1745E+01+, 
9787E+C1+, 
1175E+02+, 
3537E+02+. 
5288E+C2+. 
Of  30E  +  03+, 
9999E+99+, 


10000000EH 
86922173Eh 
70781189Ei 
55074173E. 
41117293EH 
294420U0EH 
20151870Ei 
13102586EH 
8C131869E- 
45400912E- 
23268503E- 
10368454E- 
37375932E- 
93491810E- 
10389779E- 
OOC  DOOCOE^ 


01  +  , 
00  +  , 
0C  +  , 

00  +  , 
00+, 
u0  +  , 
00+, 
00+, 

01  +  , 
01  +  , 

oi  +  , 

01  +  , 

02  +  . 

03  +  , 
03  +  , 
00+, 


oc  con 

32860 

11631 
23715 
38265 
54147 
70411 
86298 
10122 
11478 
12669 
13674 
14486 
15096 
15504 
15707 


G00E-99 
913E-01 
560E+00 
032E+00 
275E+00 
968E+00 
849E+00 
204E+0G 
996E+01 
935E+01 
033E+01 
976E+01 
008E+01 
530E+U1 
106E+01 
962E+01 
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The  sequence  of  SPS  statements  from  line  01010  to  02160 
represents  the  input  subroutine  ElljJ\/N  that  performs  a  part  of 
the  input  translations  by  converting  alphamerical  floating  point 
input  data  of  the  form  +. 50732777E-02  into  numerical  floating 
point  data  of  the  form  £073277702.   The  sequence  of  statements 
from  line  02170  to  03180  performs  an  output  translation  opposite 
to  the  above  described  input  translation.   This  output  subrou- 
tine is  labeled  EII4JN/A.   When  the  programs  A  and  B  are  to  be 
executed,  execution  starts  by  means  of  a  branch  instruction  on 
line  OI4.O2O.   If  program  switch  1  is  in  on-position,  program  A 
is  executed  according  to  the  statement  of  line  Oi+O^O.   it  is 
supposed  that  program  A  (PROG-MA)  is  executed  first.   The  5-digit 
data  field  PR0GMA+6  stores  initially  the  number  of  triples, 
(XHIGH,  MSL0PE,  BINTCP) ,  that  will  be  read-in  as  nonerasable 
information  into  core  storage;  it  also  serves  as  a  counter. 
Each  time  a  triple  (XIGH,  MSL0PE,  BINTCP)  has  been  read-in  this 
count  is  decreased  by  1.   Statement  OJ4.080  resets  this  field  to 
the  initial  value  of  l6  as  the  upper  polygonal  approximation 
consists  of  16  line  segments.   The  sequence  of  statements  from 
line  014.130  to  line  05>070  represents  the  loop  for  the  triples 
(XHIGH,  MSL0PE,  BINTCP).   At  every  pass  each  member  of  the 
triple  undergoes  the  input  translation  Ell^A/N  and  thereafter 
the  counter  PR0GMA+6  is  decremented  by  1.   The  triple  (XHIGH, 
MSL0PE,  BINTCP)  is  punched  manually  on  an  IBM  card  with  XHIGH 
starting  in  column  1  and  the  format  3E11+.8  (cf.  end  of  program). 
Statement  0f?090  causes  the  computer  to  print  out  the  message 
"ENTER  NUMBER  OP  FUNCTION  PASSES".   This  number  also  called  TIMES 
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LOOP  NUMBER  is  entered  as  a  5-digit  number  on  the  console  type- 
writer.  It  is  stored  in.  core  storage  in  field  STARTS+2.   This 
number  determines  how  many  times  ATAN(X)  has  to  be  calculated 
by  using  upper  polygonal  approximation.   The  purpose  of  repeti- 
tive calculation  of  arctan(x),  where  x  remains  unchanged,  is  to 
facilitate  computation  time  measurement.   Statement  05l50  causes 
the  console  typewriter  to  print  out  the  message  ENTER  X- .   State- 
ment 05170  causes  the  card  reader  to  read  argument  x.   Argument 
x  is  punched  manually  on  IBM  card  with  format  Eli]..  8,  e.g., 
+.11111111E+01;  and  the  mantissa  sign  is  entered  in  column  1. 
Argument  x  undergoes  likewise  the  input  translation  ElljJl/N. 
Statement  06050  transfers  the  TIMES  LOOP  NUMBER  from  field 
STARTS+6  to  count-down  area  H+6 .   The  computer  halts  at  statement 
06o5l.   Here  the  central  part  of  program  A  starts  with  the  de- 
cision tree.   This  tree  has  four  levels  and  leads  to  16  pairs 
(MSLOPE,  BINTCP) .   It  starts  at  statement  O617O  and  ends  at 
statement  120lj.O.   The  arithmetic  unit  of  this  arctangent  SPC  is 
simulated  from  12060  to  12080.   After  statement  12080  is  exe- 
cuted calculation  of  arctan(x)  is  completed  and  the  times  loop 
counter  is  decreased  by  1.   The  computer  branches  back  to  state- 
ment O6060  to  repeat  the  same  calculation.   This  looping  will 
continue  until  the  times  loop  counter  attains  zero  and  the  com- 
puter is  caused  to  halt  at  instruction  12101.   By  choosing  the 
TIMES  LOOP  NUMBER  adequately  high,  it  is  possible  to  measure 
computation  time  between  the  two  HALT  instructions  of  statement 
Olj.Olj.0  and  12101  with  a  wrist  watch.   To  cause  the  computer  to 
execute  program  B  Program  Switch  1  is  set  in  on-position.   After 
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depressing  the  Start  key  the  output  translation  Elipj/A  for 
argument  x  and  its  calculated  value  is  performed.   Both  numbers 
are  punched  with  format  2Ellj..6.   Then  the  computer  branches 
automatically  to  program  B.   The  formats  for  this  program  are 
identical  with  those  of  program  A.   The  central  part  of  this 
program  is  located  between  the  HALT  statements  131^1  and  13l8l. 
At  statement  13160  the  computer  branches  to  the  arctangent 
library  subroutine  FATN.   This  subroutine  uses  a  Chebyshev  power 
series  approximation.   Computation  time  is  measured  in  the  same 
way  as  for  program  A  and  the  same  holds  for  the  output  trans- 
lation. 

Time  measurement  yielded: 

Polygonal  approximation        28ms 
Power  series  approximation    3^4-Oms 

In  this  program  polygonal  approximation  is  twelve  times 
faster  than  power  series  approximation.   However,  this  ratio 
would  be  reduced  if  a  polygonal  approximation  is  employed  that 
yields  the  same  accuracy  of  10    as  the  power  series  approxi- 
mation. 

Higher  accuracy  for  polygonal  approximation  is  obtained  by 
employing  a  larger  number  of  line  segments.   Such  ATAN  SPC ' s 
use  more  decision  levels.   The  number  of  decision  levels  grows 
logarithmically  as  the  number  of  line  segments  increases. 
Therefore  decision  computation  time  increases  slowly,  whereas 
arithmetic  computation  time  remains  constant.   For  power  series 
approximation  computation  time  increases  exponentially  as  accuracy 
improves.   Polygonal  approximation,  therefore,  is  a  faster  com- 
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putation  method  than  power  series  approximation.   On  the  other 
hand,  polygonal  approximation  uses  more  storage  space  than 
power  series  approximation.   This  is  due  to  the  larger  number 
of  precomputed  constants  to  be  stored  and  the  greater  number 
of  instructions  involved. 
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This  thesis  discusses  a  special  purpose  digital  computer 
for  the  arctangent  function.   The  arctangent  function  is  ap- 
proximated by  upper  polygonal  approximation.   Different  kinds 
of  polygonal  approximation  are  discussed  in  detail  and  their 
features  compared.   The  "Upper  Polygonal  Approximation  with 
Uniform  Chebyshev  Norm"  turns  out  the  most  favorable  approxima- 
tion.  Problems  induced  by  different  approximations  are  thor- 
oughly treated.   Numerical  results  for  the  first  ten  approxi- 
mations are  given. 

Finally,  the  basic  ideas  for  the  realization  of  an  arc- 
tangent special  purpose  computer  are  summarized.   They  should 
convey  an  understanding  of  its  simulation  on  a  general  purpose 
computer.   Numerical  calculation  and  simulation  was  done  on  an 
IBM  1620  digital  computer. 


